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2 ■ ABSTRACT 
O '■ 

:b' 

p We present the results of a wide-field survey designed to measure the size, 

Q . inclination, and radial distributions of Kuiper Belt Objects (KBOs). The survey 

+i . found 86 KBOs in 73 square degrees observed to limiting red magnitude 23.7 

d ■ using the Canada-France-Hawaii Telescope and the 12k x 8k CCD Mosaic 

>• . camera. For the first time, both ecliptic and off-ecliptic fields were examined to 

^ . more accurately constrain the inclination distribution of the KBOs. The survey 

data were processed using an automatic moving object detection algorithm, 
allowing a careful characterization of the biases involved. In this work, we 
quantify fundamental parameters of the Classical KBOs (CKBOs), the most 
numerous objects found in our sample, using the new data and a maximum 
likelihood simulation. Deriving results from our best-fit model, we find that the 
size distribution follows a differential power law with exponent q = 4.0lo;5 (1^5 
68.27% confidence). In addition, the CKBOs inhabit a very thick disk consistent 

4-6° 

with a Gaussian distribution of inclinations with a Half- Width of ii/2 = 20°_^o 
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(1 a). We estimate that there are A^ckbos(-D > 100 km) = S.StH x 10^ (Icr) 
CKBOs larger than 100 km in diameter. We also find compelling evidence for 
an outer edge to the CKBOs at hehocentric distance R — 50 AU. 

Subject headings: Kuiper Belt, Oort Cloud — minor planets, asteroids — solar 
system: formation 

1. Introduction 

The rate of discovery of Kuiper Belt Objects (KBOs) has increased dramatically since 
the first member (1992 QBi) was found (Jewitt & Luu 1993). As of Dec 2000, ~ 400 KBOs 
are known. These bodies exist in three dynamical classes (Jewitt, Luu & Trujillo 1998): 
(1) the Classical KBOs (CKBOs) occupy nearly circular (eccentricities e < 0.25) orbits 
with semimajor axes 41 AU ^ a ^ 46 AU, and they constitute ~ 70 % of the observed 
population; (2) the Resonant KBOs occupy mean-motion resonances with Neptune, such as 
the 3:2 (the Plutinos, a ^ 39.4 AU) and 2:1 (a ^ 47.8 AU), and comprise ~ 20 % of the 
known objects; (3) the Scattered KBOs represent only ~ 10 % of the known KBOs, but 
possess the most extreme orbits, with median semimajor axis a ~ 90 AU and eccentricity 
e ~ 0.6, presumably due to a weak interaction with Neptune (Duncan & Levison 1997, 
Luu et al. 1997, and Trujillo, Jewitt & Luu 2000). Although these classes arc now well 
established, only rudimentary information has been collected about their populations. One 
reason is that only a fraction of the known KBOs were discovered in well-parametrized 
surveys that have been published in the open literature (principally Jewitt & Luu 1993 (1 
KBO); Jewitt & Luu 1995 (17 KBOs); Irwin, Tremaine & Zytkow 1995 (2 KBOs); Jewitt, 
Luu & Chen 1996 (15 KBOs); Gladman ct al. 1998 (5 KBOs); Jewitt, Luu & Trujillo 
1998 (13 KBOs); and Chiang & Brown 1999 (2 KBOs)). In this work, we characterize 
the fundamental parameters of the CKBOs: the size distribution, inclination distribution 
and radial distribution using a large sample (86 KBOs) discovered in a well-characterized 
survey. 

The quintessential measurement of the size distribution relies on the Cumulative 
Luminosity Function (CLF). The CLE describes the number of KBOs deg^^ (S) near 
the ecliptic as a function of apparent red magnitude {rriji). It is fitted with the relation 
log(S) = a{mji — mo), where mo is the red magnitude at which E = 1 KBO deg~^. The 
slope (a) is related to the size distribution (described later). Although many different 
works have considered the CLF, two papers are responsible for discovering the majority of 
KBOs found in pubhshed surveys: Jewitt, Luu & Chen (1996) and Jewitt, Luu & Trujillo 
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(1998). The former constrained the CLF over a 1.6 magnitude range (23.2 < mu < 24.8) 
with 15 discovered KBOs while the latter covered a complimentary 2.5 magnitude range 
(20.5 < rriji < 23.0), discovering 13 objects. Jewitt, Luu & Trujillo (1998) measured the 
CLF produced from these two data sets and found a = 0.58 ± 0.05 and mo = 23.27 ± 0.11. 
Gladman et al. (1998) criticized this work on 2 main counts: (1) they believed that Jewitt, 
Luu & Chen (1996) underestimated the number of KBOs and (2) the fit in the Jewitt, Luu 
& TrujiUo (1998) survey used a least-squares approach that assumed Gaussian errors rather 
than Poissonian errors. Gladman et al. (1998) found 5 additional KBOs and re-analysed 
the CLF using a Poissonian maximum hkelihood method to refit the CLF to (1) the Jewitt, 
Luu & Trujillo (1998) data without the Jewitt, Luu & Chen (1996) data and (2) a fit to 
the 6 different surveys available at the time except for Tombaugh (1961), Kowal (1989) and 
Jewitt, Luu & Chen (1996). Both of these fits were steeper but formally consistent with 
the original Jewitt, Luu & Trujillo (1998) data at the ~ 1.5a level: (1) a = 0.72liQ'^Q and 
mo = 23.3l[]i and (2) a = 0.76l|]:}? and mo = 23AOtoil- Chiang & Brown (1999) find a 
flatter size distribution of a = 0.52 ± 0.05 and mo = 23.5 ± 0.06 much closer to the Jewitt, 
Luu & Trujillo (1998) result. They observed that the steep size distribution reported by 
Gladman et al. (1998) was an artifact of their selective exclusion of part of the available 
survey data, not of their use of a different fitting method. The first goal of this work is to 
measure the CLF and additionally constrain the power-law slope of the size distribution 
using a single well-characterized survey and a maximum likelihood simulation which allows 
for the correction of observational biases. 

An accurate characterization of the inclination distribution of the KBOs is critical to 
understanding the dynamical history of the outer Solar System since the era of planetesimal 
formation. We expect that the KBOs formed by accretion in a very thin disk of particles 
with a small internal velocity dispersion (e.g. Kenyon & Luu 1998 and Hahn & Malhotra 
1999) and a correspondingly small inclination distribution. However, the velocity dispersion 
indicated by the inclination distribution in the present-day Kuiper Belt is large. Jewitt, 
Luu & Chen (1996) measured the apparent Half- Width of the Kuiper Belt inclination 
distribution to be ~ 5°. They noted a strong bias against observing high inchnation objects 
in ecliptic surveys, and they estimated the true distribution to be much thicker, with an 
inclination distribution Half- Width of ^ 15°, corresponding to a vertical velocity dispersion 
of ~ 1 km/s. Several conjectures have been advanced to explain the thickness of the Kuiper 
Belt: Earth-mass planetesimals may have been scattered through the belt in the late stages 
of the planet-formation era, exciting the Kuiper Belt (Morbidelli & Valsecchi 1997 and Petit 
et al. 1999); stellar encounters may have enhanced the velocity dispersion of the distant 
KBOs (Ida, Larwood & Burkert, 2000); and the dispersion velocity of small bodies tends to 
grow to roughly equal the escape speed of the bodies contributing the most mass (the large 
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bodies for size distributions with < 4) in the belt (Aarseth, Lin & Palmer 1993). As there 
is much speculation about the origin of the large velocity dispersion of the Kuiper Belt, but 
only one published measurement (Jewitt, Luu & Chen 1996), the second goal of this work 
is to accurately quantify the inclination distribution from our large sample of objects. 

The radial extent of the Classical Kuiper Belt has not been well constrained. None of 
the CKBOs have been discovered beyond 50 AU. This trend was first noted by Dones 
(1997) who suggested that the 50-75 AU region may be depleted; he found the results of 
a Monte- Carlo simulation of CKBOs drawn from a rather flat differential size distribution 
(power-law index g = 3) to be inconsistent with the observations of the 6 CKBOs discovered 
by Jewitt, Luu & Chen (1996). Jewitt, Luu & Trujillo (1998) discovered all of their KBOs 
at heliocentric distances i? < 46 AU. In the absence of other effects, one should expect 
to find fewer bodies with i? > 50 AU than with i? ~ 40 AU, as the former are about a 
magnitude fainter than the latter. However, through the use of a Monte-Carlo model they 
demonstrated that the bias against objects beyond 50 AU is not strong enough to explain 
the distribution of discovery distances. They speculated that the lack of bodies discovered 
beyond 50 AU could be caused by a combination of (1) a decrease in the maximum KBO 
size (and reduction in the brightest and most detectable objects) beyond 50 AU or (2) the 
size distribution might steepen beyond 50 AU, putting more of the mass in the smaller, 
less-detectable bodies. They also suggested that the lack of i? > 50 AU objects could be 
explained by an outer edge to the Classical Kuiper Belt at 50 AU. 

Two later papers questioned the existence of an edge to the Kuiper Belt near 50 AU. 
Gladman et al. (1998) suggested that the number of objects expected to be discovered 
beyond 50 AU is highly dependent on the size distribution because steep size distributions 
reduce the number of large (bright) bodies relative to small (faint) bodies. Gladman et 
al. (1998) adopted a relatively steep distribution (q — 4.65), and found no significant 
evidence of a truncated belt. Chiang &: Brown (1999) found that 8%-13% of the ~ 100 
objects known at the time should have been found beyond 50 AU, and suggested that this 
precludes the presence of a density enhancement beyond 50 AU, but could not definitively 
rule out a density deficit. Allen, Bernstein & Malhotra (2001) have also recently reported 
the detection of an outer edge to the Kuiper Belt. The third goal of the present work is 
to test the distribution of the discovery distances for the presence of an outer edge to the 
Kuiper Belt. 
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2. Survey Data 

Observations were made at the 3.6 m diameter Canada-France-Hawaii Telescope using 
the 12288 x 8192, 15 /xm pixel mosaic CCD (CFHT 12k; Cuillandre et al. 2000). Built 
at the University of Hawaii (UH) the CFHT 12k comprises 12 edge-abutted, thinned, 
high quantum efficiency (QE ~ 0.75), 4096 x 2048 pixel Lincoln Laboratory CCDs. It is 
currently the largest close-packed CCD camera in the world. When mounted at the CFHT 
f/4 prime focus the camera yields a plate scale of 0.206 arc sec/pixel, corresponding to a 
0.330 sq deg field of view in each 200 Mbyte image. Images were taken through a Mould 
R filter, with a central wavelength of 6581 A and a bandwidth of 1251 A. Instrumental 
parameters of the survey are summarized in Table |^. 

Observations were taken within a few days of new moon under photometric conditions 
during three periods: Feb 10 - 15 1999, Sep 5-8 1999, and Mar 31 - Apr 3 2000. Fields 
were imaged at airmasses < 1.7 and were within 1.5 hours of opposition. We chose to use 
short 180 s exposures at the CFHT to maximize area coverage and detection statistics. All 
discovered objects were accessible for recovery at the UH 2.2 m telescope during comparable 
seeing conditions with exposure times of < 600 seconds. Each field was imaged three times 
(a "field triplet"), with about 1 hour timebase between exposures. Fields imaged appear in 
Figures |I], |^ and ^ and in Table ^. The CFHT observations were taken at three ecliptic 
latitudes (3 = 0°, 10°, and 20° to probe the inclination distribution of the KBOs (see §|p. 

Photometric calibrations were obtained from Landolt (1992) standard stars imaged 
several times on each chip. Three CFHT 12k chips of poor quality were replaced between 
the Feb 1999 and Sep 1999 runs. The positions of four other CFHT 12k chips within the 
focal plane array were changed to move the cosmetically superior chips towards the center 
of the camera. The photometric calibration accounts for these changes, as shown in Table ^ 
containing the measured photometric zero points of the chips. In addition, chip 6 was not 
used in Feb 1999 because of its extremely poor cosmetic quality. The area covered in the 
fields from Feb 1999 was corrected for this 8% reduction in field-of-view. The area imaged 
in Mar 2000 included some small field overlap (6%), resulting in a minor correction applied 
to the reported total area surveyed. 

Each of the 12 CCDs in the CFHT 12k functions as an individual detector, with 
its own characteristic bias level, flat field, gain level, and orientation (at the ~ 1° level). 
The bias level for each chip was estimated using the row-by-row median of the overscan 
region. Flat fields were constructed from a combination of (1) the median of normalized 
bias-subtracted twilight fiat fields and (2) a median of bias-subtracted data frames, with a 
clipping algorithm used to remove excess counts due to bright stars. Fields were analysed 
by subtracting the overscan region, dividing by the composite fiats and searching for moving 
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objects using our Moving Object Detection Software (MODS, Trujillo & Jewitt 1998). We 
rejected bad pixels through the use of a bad pixel mask. 

Artificial moving objects were added to the data to quantify the sensitivity of the 
moving object detection procedure (Trujillo & Jewitt 1998). The seeing during the survey 
typically varied from 0.7 arc sec to 1.1 arc sec (FWHM). Accordingly, we subdivided and 
analysed the data in 3 groups based on the seeing. Artificial moving objects were added to 
bias-subtracted twilight sky-flattened images, with profiles matched to the characteristic 
point-spread function for each image group. These images were then passed through the 
data analysis pipeline. The detection efficiency was found to be uniform with respect to 
sky-plane speed in the 1-10 arc sec/hr range. At opposition, the apparent speed in arc 
sec/hr, 6*, of an object is dominated by the parallactic motion of the Earth, and follows 
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where R is heliocentric distance in AU (Luu & Jewitt 1988). From Equation |1], our 
speed limit criterion for the survey, 1 " hr^^ < ^ < 10 " hr~^, corresponds to opposition 
heliocentric distances 10 AU <^ R <^ 140 AU, with efficiency variations within this range 
due only to object brightness and seeing. 



ftanh r"^°-"^^ )+l), (2) 



The magnitude-dependent efficiency function was fitted by 

( . 

where < £ < 1 is the efficiency with which objects of red magnitude m^j are detected, emax 
is the maximum efficiency, m^jso is the magnitude at which e = emax/2, and a magnitudes 
is the characteristic range over which the efficiency drops from ^max to zero. Table § shows 
the efficiency function derived for each seeing category, along with an average of the seeing 
cases, weighted by sky area imaged, applicable to the entire data set. The efficiency function 
is known to greater precision than the ~ 0.1 magnitude uncertainty on our discovery 
photometry. Changes to the efficiency function of < 0.1 magnitudes produce no significant 
variation in our results for the size or inclination distributions. 

The MODS software, running on two Ultra 10 computers, was fast enough to efficiently 
search for the KBOs in near real-time, so that newly detected objects could be quickly 
discovered and re-imaged. We imaged ~ 35 field triplets each night at the CFHT, 
corresponding to ~ 20 Gbytes of raw data collected per night, plus several more Gbytes 
for flat flelds and standard stars. Eighty-six KBOs were found in the CFHT survey, 2 of 
which were serendipitous re-detections of known objects. The discovery conditions of the 
detected objects appear in Table H. Photometry was performed using a 2.5 arc second 
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diameter synthetic aperture for discovery data, resulting in median photometric error of 
0.15 magnitudes and a maximum photometric error of 0.3 for the faintest objects. Our 
results are unaffected by this error; randomly introducing ±0.15 magnitude errors in our 
simulations (described later), and ±0.3 magnitude errors in the faintest objects produced 
no statistically significant change. Trailing loss was insignificant as the KBOs moved only 
0.15 arc sec during our integration. 

2.1. Recovery Observations and Orbits 

Extensive efforts were made to recover all objects using the UH 2.2 m telescope. 
Attempts were made to recover the objects one week after discovery, then one, two and 
three months after discovery. Most of these attempts were successful, as demonstrated 
by the fact that 79 of the 86 CFHT objects were recovered. The loss of 7 objects is the 
result of unusually poor weather during the Mar-May 1999 recovery period. Only 6 of 
the 79 recovered objects have arc-lengths shorter than 30 days as of Dec 1, 2000. Second 
opposition observations have been acquired for 36 of the 78 KBOs found in Feb 1999 and 
Sep 1999. 

Orbits derived from the discovery and recovery data appear in Table ^. The listed 
elements are those computed by Brian Marsden of the Minor Planet Center. We also 
benefited from orbital element calculations by David Tholen (Univ. of Hawaii). Both 
sources produced comparable orbital solutions to the astrometric data. 

With only first opposition observations, the inclination and heliocentric distance at 
discovery can be well determined for nearly all KBOs, as depicted in Figure ^. We find that 
the semimajor axis and eccentricity determinations are less reliable but are usually good 
enough to classify the objects as either Classical, Resonance or Scattered KBOs, as depicted 
in Figure |. We find that 6 out of 36 (17%) of the objects exhibit orbital changes large 
enough for their dynamical classification to change from the first opposition to the second 
opposition. Randomly rejecting 17% of our sample (to simulate misclassification) does not 
significantly change the results. In addition, rejection of all but the multi-opposition objects 
does not significantly change our results; as expected, the total number of KBOs estimated 
decreased by a factor ~ 2 and error bars increased by a factor ~ \pi due to the sample size 
reduction. The eccentricity and semimajor axes of all objects with a < 50 AU (this includes 
all Classical KBOs) appear in Figure 

In the next two sections, we use our observations to constrain three fundamental 
quantities of the Classical KBOs: (1) the size distribution index, g, (2) the half-width of 
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the inclination distribution, and (3) the total number of CKBOs larger than 100 km 
in diameter, A^^ckbosI-D > 100 km). The quantities ii/2 and q are uncorrected, as the 
observable constraining ii/2 is the inclination distribution and the observable constraining 
q is the absolute magnitude distribution. However, iVcKBOs(-D > 100 km) is a function of 
both q and ii/2 as a steeper size distribution or thicker inclination distribution will each 
allow more bodies to be present. In the maximum likelihood simulations that follow, the 
ideal case would be to constrain q, ii/2 and iVcKBOs(-D > 100 km) and estimate errors in one 
simulation, however, this is difficult computationally. Therefore, we find the best-fit values 
of the three parameters in a single simulation, but estimate the errors on the parameters in 
two simulations, one that estimates the q-NcKBOsiD > 100 km) joint errors and one that 
estimates the «i/2-A'^ckbos(-D > 100 km) joint errors. We then combine the two simulation 
results in quadrature to determine the errors on A''ckbos(-D > 100 km). 

3. The Size Distribution of the Classical KBOs 

We estimate the size distribution of the KBOs from our data in two ways. The first is 
a simple estimate made directly from the distribution of ecliptic KBO apparent magnitudes 
(CLF). The second is a model which simulates the discovery characteristics of our survey 
through the use of a maximum likelihood model constrained by the absolute magnitude of 
the Classical KBOs. 



3.1. Cumulative Luminosity Function 

We model the CLF with a power-law relation, logS = a(m^ — mo) (§|T])- The KBOs 
are assumed to follow a differential power-law size distribution of the form n{r)dr oc r~'^dr, 
where n{r)dr is the number of objects having radii between r and r + dr, and q is the index 
of the size distribution. Assuming albedo and heliocentric distance distributions that are 
independent of KBO size, the simple transformation between the slope of the CLF (a) and 
the exponent of the size distribution (g) is given by (Irwin et al. 1995), 

g = 5a + l. (3) 

Under these assumptions, the size distribution can be estimated directly from the CLF. 

We estimated the CLF by multiplying the detection statistics from the observed 
distribution of object brightnesses by the inverse of the detection efficiency. We assumed 
Poisson detection statistics, with error bars indicating the interval over which the integrated 
Poisson probability distribution for the observed number of objects contains 68.27% of 
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the total probability (identical to the errors derived by Kraft, Burrows & Nousek 1991). 
This is nearly equal to the Gaussian case for all data points resulting from more than a 
few detections. We have included all 74 KBOs discovered in our 37.2 sq deg of ecliptic 
fields in the estimate of the CLF. This includes the lost objects, as the CLF is simply 
a count of the number of bodies discovered at a given apparent magnitude. Our results 
appear in Figure 0, with other published KBO surveys. All observations were converted 
to -R-band if necessary assuming V — R = 0.5 for KBOs (Luu & Jewitt 1996), and error 
bars were computed assuming Poisson detection statistics. The data point of Cochran et 
al. (1995) near m/j = 28 was omitted because of major uncertainties about its reliability 
(Brown, Kulkarni & Liggett 1997, cf. Cochran et al. 1998). Early photographic plate 
surveys (Tombaugh 1961, portions of Luu & Jewitt 1988, and Kowal 1989) have unproven 
reliability at detecting faint slow-moving objects, and plate emulsion variations and defects 
make accurate photometric calibration difficult. The photographic plate survey data were 
not used in our analysis. 

The CLF points are highly correlated with one another, resulting in a heavy weighting 
of the bright object data points. Thus, we fitted the Differential Luminosity Function 
(DLF) instead. We plot the DLF data points at the faint end of the bin, representing 
the modal value in that bin. Very small bin sizes were chosen (0.1 magnitudes) to negate 
binning effects incurred from averaging the detection efficiency (Equation ^ over a large 
magnitude range. For any non-zero CLF slope, a, the DLF and CLF slopes are equal due 
to the exponential nature of the CLF. The DLF was modelled by evaluating the Poisson 
probability of detecting the observed DLF given a range of mo and a, with the maximum 
probability corresponding to our best-fit values. Error bars were determined by finding 
the contours of constant joint probability for mo and a enclosing 68.27% of the total 
probability, a procedure similar to that used below for the maximum likelihood simulation. 
Computations from this procedure are summarized in Table |^. We find that the slope 
of the CLF is a = 0.64l|]:}J with mo = 23.23l[ji^, which corresponds to g = 4.2 ± 0.5 
from Equation |^. We also fitted the CLF by applying the maximum likelihood method 
described by Gladman et al. (1998) to our data, which yields statistically identical results 
to the binned DLF procedure: a = 0.63 ± 0.06 and mo = 23.04lo;og, corresponding to 
q = 4.2 ± 0.3. The maximum-likelihood method provides slightly better signal-to-noise, and 
is independent of binning effects. However, it does not provide a visualization of the data, 
as presented for the DLF fit. We adopt the maximum-likelihood procedure as our formal 
estimate of the CLF slope. Both methods estimating the size distribution are in statistical 
agreement with the more detailed analysis presented in the next section. 

The best-fit a = 0.63 magnitude distribution was compared with the observed 
magnitude distribution using a Kolmogorov-Smirnov test (Press et al. 1992), producing 
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a value of D = 0.13. If the model and the data distributions were identical, a deviation 
greater than this would occur by chance 12% of the time. Thus, our linear model is not a 
perfect fit, but it is statistically acceptable. 

3.2. Maximum Likelihood Simulation 

We now present more detailed analysis of the size distribution. Since we model the 
detection statistics of an assumed population, we choose to model the 49 Classical KBOs 
discovered on the ecliptic as they are numerically dominant in the observations and their 
orbital parameters are more easily modelled than other KBO classes. Our selection criteria 
for CKBOs are perihelion q' > 37 AU and 40.5 AU < a < 46 AU. Given the size of an 
object and its orbital parameters, we can compute its position, velocity, and brightness, 
allowing a full "Monte-Carlo" style analysis of the bias effects of our data collection 
procedures. The apparent brightness was computed from: 

m = mQ- 2.5 log(p<l>(a>2) + 2.5 log(2.25 x lO^^i^^A^), (4) 

where a' is the phase angle of the object, is the Bowell et al. (1989) phase function, 

geometric albedo is given by p, r is the object radius in kilometers, R is the heliocentric 
distance, and A is the geocentric distance, both in AU (Jewitt & Luu 1995). The apparent 
red magnitude of the Sun was taken to be = —27.1. For this work, we assume p = 0.04, 
consistent with a Centaur-like albedo (Jewitt & Luu 2000). We neglect phase effects 
(setting $(«') = 1) since the maximum phase angle of an object at i? = 40 AU within 1.5 
hours of opposition is a' = 0.55°. This corresponds to $(«') = 0.91, a change in brightness 
of only 0.09 magnitudes, which is less than other uncertainties in the data. 

This apparent brightness is used in a biasing-correction procedure (Trujillo, Jewitt & 
Luu 2000 and Trujillo 2000), summarized here: 

1. A model distribution of KBOs is assumed (described in Table H). 

2. KBOs are drawn randomly from the model distribution. 

3. For each KBO, the apparent speed and ecliptic coordinates are computed from the 
equations of Sykes Sz Moynihan (1996, a sign error was found in Equation 2 of their 
text and corrected), and compared to the observed fields and speed criteria. 

4. The apparent magnitude is computed from Equation ^. 

5. The efficiency function (Equation ^ and our field area covered are used to determine 
if the simulated object would be "detected" in our survey. 
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6. A histogram of the detection statistics for the simulated objects is constructed, 
logarithmically binned by object size for the size distribution model and binned by 
inclination for the inclination-distribution model. Binning effects were negligible due 
to small bin choice. 

7. Steps 1-6 are repeated until the number of detected simulated objects is at least a 
factor 10 greater than the number of observed objects in each histogram bin (typically 
requiring a sample of 10® < N < 10® simulated objects, depending on the observed 
distribution). 

8. The likelihood of producing the observed population from the model is estimated by 
assuming that Poisson detection statistics (P = ^ exp (— /i)) apply to each histogram 
bin, where /i represents the expected number of simulated objects "discovered" given 
the number of objects simulated and n represents the true number of KBOs observed. 
Thus, the observed size distribution, calculated from Equation ^ is used to constrain 
the q model, and the observed inclination distribution is used to constrain the i model 



These steps are repeated for each set of model parameters in order to estimate the likelihood 
of producing the observations for a variety of models. 

For the size distribution analysis, we take our best-fit model of the width of the 
inclination distribution (Half- Width ii/2 = 20°, as estimated in the next section), and vary 
the size distribution index q, and the total number of objects A^ckbos(-D > 100 km). Model 
parameters are summarized in Table || and results appear in Figure ||. Our best-fit values 
are 

q = A.Otoi (1^) and 

= 4.0t^:? (3a), 
and 

iVcKBOs(P' > 100 km) = 3.8tli x 10^ (la) and 

= 3.8^1^ X 10^ (3a), 

where the errors for A'"ckbos(-D > 100 km) have been combined in quadrature from the 
results of the q and ii/2 fits, as described at the end of § |2]T|. The values for q are consistent 



with previously published works (Table |^) and the q derived from the CLF data in the 
simple model (Equation |^). The results are consistent with the distribution of large 
(-D > 150 km) main-belt asteroids {q = 4.0, Cellino, Zappala, & Farinella 1991) and rock 
crushed by hypervelocity impacts (g = 3.4, Dohnanyi 1969). In addition, the scenario 
where the cross-sectional area (and thus optical scattered light and thermal emission) 
is concentrated in the largest objects (g < 3, Dohnanyi 1969) is ruled out at the > 2a 
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(> 95.4% confidence) level. Our results are also consistent with Kenyon & Luu (1999) who 
simulate the growth and velocity evolution of the Kuiper Belt during the formation era in 
the Solar System. They find several plausible models for the resulting size distribution, all 
of which have g ^ 4. In Figure ^ we plot the best-fit model CKBO distribution with the 
observed DLF to demonstrate the expected results from different size distributions. 

The magnitude distribution expected from the maximum likelihood model was 
compared to the observed magnitude distribution, as was done for the CLF-derived 
magnitude distribution in § p.l| . The Kolmogorov-Smirnov test produced D = 0.17; a 
greater deviation would occur by chance 11% of the time. 

In our Classical KBO maximum likelihood simulation, we have ignored possible 
contributions of the 7 lost KBOs, since their orbital classes are not known. However, 
including them in the simulations by assuming circular orbits at the heliocentric distance 
of discovery results in statistically identical resulted for q, and the expected 7/49 rise in 
A^CKBOs(/^ > 100 km). 



4. Inclination Distribution of the Classical KBOs 

The dynamical excitation of the Kuiper Belt is directly related to the inclination 
distribution of the KBOs. We present the inclinations of the CKBOs found in the CFHT 
survey in Figure [ly. Assuming heliocentric observations, a KBO in circular orbit follows 

sin f3 = sin i sin / (5) 

where /? is the heliocentric ecliptic latitude, < i < 90° is the inclination, and < / < 360° 
represents the true anomaly of the object's orbit with / = and 180° representing the 
ecliptic plane crossing (the longitude of perihelion is defined as in this case). Using 
Equation ||, we plot the fraction of each orbit spent at various ecliptic latitudes as a function 
of i (Figure |lT]). This plot demonstrates two trends concerning the ecliptic latitude of 
observations /?obs- First, high-inclination objects are a factor 3-4 times more likely to be 
discovered when /Sobs ~ i than when observing at low echptic latitudes {Pohs < i)- Second, 
the number of expected high- inclination objects drops precipitously, roughly as once 
i > 1.5/3obs (Jewitt, Luu & Chen 1996). 

These facts led us to observe at three different ecliptic latitudes (0°, 10° and 20°) 
to better sample the high-inclination objects. During two observation periods (Sep 1999 
and Mar 2000) care was made to interleave the ecliptic fields with the off-ecliptic fields 
on timescales of ~ 30 minutes. This technique provides immunity to drift in the limiting 
magnitude which might otherwise occur in response to typical slow changes in the seeing 
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through the night. The results for the robust, interleaved fields matched those for the 
seeing-corrected Feb 1999 fields where fields were interleaved on much longer timescales 
of ~ 3 hours. Accordingly, we combined the data sets from all epochs to improve 
signal-to-noise. In the next sections, we analyse the inclination distribution using two 
techniques to demonstrate the robustness of our method. 



4.1. Simple Inclination Model 

First, since fields were imaged at three different ecliptic latitudes, the surface density of 
objects at each latitude band (S(0°), S(10°) and S(20°)) can directly yield the underlying 
inclination distribution. In our simple model, we generate an ensemble of inclined, 
circular orbits drawn from a Gaussian distribution centered on the ecliptic, and having 
a characteristic Half- Width of zi/2. The probability of drawing a KBO with inclination 
between i and i + di is given by 



where a = «i/2/V2 ln2. Using this relation, and Equation |^, we simulate the expected 
values of S(0°), S(10°) and S(20°) for various ii/2- These are compared to two ratios 
measured from our observations, i?(10°,0°) = S(10°)/S(0°) and i?(20°,0°) = S(20°)/S(0°) 



Results appear in Table |T0|, and demonstrate that the characteristic half-width of the 



inclination distribution in the Kuiper Belt is ii/2 ~ 17°~^^o (Icr = 68.27% confidence). This 
simple model does not use the observed inclination distribution of the individual objects, 
merely the surface density of objects found at each ecliptic latitude, thus we have combined 
all objects from all KBO classes into this estimate. 



4.2. Full Maximum Likelihood Inclination Model 



Second, we use the maximum likelihood model described in 



We list the 



parameters of the model in Table |Tl|. This model encompasses the additional constraint 
of the observed inclination distribution, as well as the parallactic motion of the Earth and 



KBO orbital motion to produce more realistic results. Results appear in Figure |12| with 
iVcKBOs(-D > 100 km) representing the number of CKBOs with diameters greater than 100 
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km. The maximum likelihood occurs at 



2^/2 = 20°+5o (1(t) and 



^^-10+26 

and 



20°+p (3or), 



iVcKBOs(^ > 100 km) = 3.8tli x 10^ (la) and 

= 3.8l|^ X 10^ (3a), 

where the errors for Ai'cKBOs(-D > 100 km) have been estimated from the ii/2 and q fits, 
combined in quadrature, as described at the end of § |2.1| . This maximum likelihood model 
is consistent with the simple model described in § [4.1| . In Figure |TB|, we plot the observed 



surface density of objects as a function of ecliptic latitude and compare these data to our 
best-fit models. This illustrates the fundamental fact that even though the true inclination 
distribution of the KBOs is very thick (ii/2 ~ 20°), the surface density drops off quickly 
with ecliptic latitude, reaching half the ecliptic value at an ecliptic latitude of /3 ~ 3° 
(E(3°)/S(0°) < 0.5). 

The functional form of the inclination distribution cannot be well constrained by our 
data. However, the best-fit Gaussian distribution was compared to a fiat-top ("top-hat") 
inclination distribution, with a uniform number of objects in the 0° < i < 30° range. The 
Gaussian and fiat-top models were equally likely to produce the observed distribution in 
the 65% confidence limit (< la). A Gaussian model multiplied by sin(z) was also tried but 
could be rejected at the > 3a level because it produced too few low-inclination objects. 
We also tested the best-fit model presented by Brown (2001), consisting of two Gaussians 
multiplied by sin(z), 

[aexp(^) + (l-a)exp(^)]sinz, (7) 

where a = 0.93, ai = 2.2°, and a2 = 18°, and found it equally compatible with our single 
Gaussian model (Equation P). Because the Gaussian model was the simplest model that 
fit the observed data well, we chose it to derive the following velocity dispersion results. 

We first find the mean velocity vector of all the simulated best-fit CKBOs, v, in 
cylindrical coordinates (normal vectors f, 9, and z representing the radial, longitudinal 
and vertical components respectively). The mean velocity vector v is consistent with a 
simple Keplerian rotation model at i? ~ 46 AU. We then compute the relative velocity of 
each KBO from this via It; — ■i'il, where vl is the velocity dispersion contribution of the ith 
KBO. We find the resulting root-mean-square (RMS) velocity dispersion of the f, 6, and 
z components to be equal to Avr = 0.51 km/s, Avg = 0.50 km/s, and Av^ = 0.91 km/s, 
combining in quadrature for a total velocity dispersion of Av = J Av^ + Avg + AvI = 1.16 
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km/s. An error estimate of the velocity dispersion can be found by following a similar 
procedure for the = 16°and 26°(±lcr) models, yielding Av = LlGt^ie km/s. 



4.3. Inferred Mass 

The Kuiper Belt mass inferred from these results can be directly calculated from the 
size distribution and the number of bodies present. For the best-fit g = 4.0 size distribution, 
the mass of CKBOs in bodies with diameters -Dmin < D < -Dmax is 

4 

max / -7rprin(L'rnax/Anm), (8) 

where p is the bulk density of the object. The normalization constant F is calculated from 
the results of our simulation, 

F ^ 3.0 X 10^2 Yii^p-^^-^N{D > 100 km), (9) 

where N{D > 100 km) = 3.8 x 10^ (§13), yielding F = 1.4 x 10^^ m^ assuming pr = 0.04. 
The mass for 100 km < D < 2000 km then becomes 

M(100 km, 2000 km) ^ 0.03Mm ( ^ ) ( 

^ ' ' ^ V 1000 kg m-^ J \ 

where = 6.0 x 10^^ kg is the mass of the earth. The uncertainties on this value are 
considerable as the characteristic albedo and density of the CKBOs are unknown. 




4.4. Comparison of the Classical KBOs to Other Dynamical Classes 

We found that the total number of CKBOs is given by A/'ckbosI-D > 100 km) = 
3.8li;5 ^ 10^- This can be compared to the other main dynamical populations (the Resonant 
and Scattered KBOs) from our data. Observational biases favor the detection of the 
Plutinos over the Classical KBOs due to their closer perihelion distance. We found only 7 
Plutinos (4 ecliptic and 3 off-ecliptic) so we can make only crude (factor ~ 2 statements) 
about the true size of the population. Thus, we use the results of Jewitt, Luu & Trujillo 
(1998) who estimate that the apparent fraction of Plutinos {Pa) in the Kuiper Belt is 
enhanced relative to the intrinsic fraction (Pj) by a factor Pa/ Pi ~ 2 for g = 4.0 and 
^max = 1000 km. Applying this correction to our ecliptic observations (4 Plutinos and 49 
Classical KBOs) indicates that the total number of Plutinos larger than 100 km in diameter 
is quite small, 

N^iutuUD > 100 km) ^ -^:^iVcKBOs ~ MOO. (11) 



- 16 - 



The populations of the Plutinos and the 2:1 Resonant objects are important measures 
of the resonance sweeping hypothesis (Malhotra 1995), which predicts equal numbers of 
objects in each resonance. Since the 2:1 objects are systematically farther from the sun 
than the Plutinos, the true Plutino/2:l ratio is higher than the observed ratio. Jewitt, 
Luu & Trujillo (2000) estimate the observed/true bias correction factor to be ~ 0.310 for a 
survey similar to ours (g = 4 and rriR^Q = 24.0). Only 2 of our objects (both found on the 
ecliptic) are < 0.5 AU from the 2:1 Resonance, so we find the Plutino/2:l fraction is given 
by (4/2)0.310 ~ 0.6. Due to the small number of bodies involved, this is only an order 
of magnitude estimate. Within the uncertainties, our observations are consistent with the 
hypothesis that the 3:2 and 2:1 resonances are equally populated. 

The observational biases against the Scattered KBOs are considerable. Trujillo, 
Jewitt & Luu (2000) estimate the total population of the Scattered KBOs to be 
A'sKBOs(-D > 100 km) = 3.ltJ;3 x 10"^, approximately equal to the population of Classical 
KBOs derived from our data. We summarize the relative populations by presenting their 
number ratios: 

Classical : Scattered : Plutmo : Resonant 2:1 = 1.0 : 0.8 : 0.04 : 0.07. (12) 



5. The Edge of the Classical Kuiper Belt 

We found no objects beyond heliocentric distance i?obs = 48.9 AU. There are two 
possibilities to explain this observation: (1) this is an observational bias effect and the 
bodies beyond i?obs cannot be detected in our survey, or (2) there is a real change in the 
physical or dynamical properties of the KBOs beyond -Robs- In order to test these two 
possibilities, we compare the expected discovery distance of an untruncated Classical Kuiper 



Belt to the observations, as depicted in Figure III. This untruncated CKBO distribution 



is identical to our best-fit model from § [4.2|, except 40.5 AU < a < 200 AU, instead 



of 40.5 AU < a < 46 AU. The total number of bodies produced was considered a free 



parameter in this model. Inspecting Figure 0, the absence of detections beyond 50 AU is 
inconsistent with an untruncated model with radial power to the ecliptic plane surface 
density. Assuming Poisson statistics apply to our null detection beyond -Rmax, the 99.73% 
(3(j) upper limit to the number of bodies (/i) expected beyond -Rmax can be calculated 
from 1 — 0.9973 = exp(— /x), yielding = 5.9 KBOs. We found 49 ecliptic Classical KBOs 
inside the -Rmax limit, so the 3a upper limit to the number density of KBOs beyond -Rmax 
is 49/5.9 ~ 8 times less than the number density of Classical KBOs. Although we have 
constrained the outer edge by the heliocentric distance at discovery -R, which is a directly 
observable quantity, a dynamical edge would be set by the semimajor axes (a) of the object 
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orbits. This difference has httle effect on our findings as the known CKBOs occupy nearly 
circular orbits with median eccentricity e = 0.08 (the calculated median is conservative 
as it includes only bodies with e > to protect against short-arc orbits, which typically 
assume e = 0). Since an untruncated distribution (1) is incompatible with our data, we 
must conclude that scenario (2) applies — there must be a physical or dynamical change in 
the KBOs beyond -Rmax- 

There are several possible physical and dynamical mechanisms that could produce the 
observed truncation of the belt beyond -Rmax = 50 AU (Jewitt, Luu & Trujillo 1998): (1) 
the size distribution of the belt might become much steeper beyond -Rmax, putting most of 
the mass of the belt in the smallest, undetectable objects; (2) the size distribution could be 
the same (g = 4), but there might be a dearth of large (i.e. bright) objects beyond -Rmax, 
suggesting prematurely arrested growth; (3) the objects beyond -Rmax may be much darker 
and therefore remain undetected; (4) the eccentricity distribution could be lower in the 
outer belt, resulting in the detection of fewer bodies; (5) the ecliptic plane surface density 
variation with radial distance may be steeper than our assumed p = 2; and (6) there is a 
real drop in the number density of objects beyond -Rmax- We consider each of these cases in 
turn, and their possible causes. 

Detailed simulations of the growth of planetesimals in the outer Solar System have 
not estimated the radial dependence of the formation timescale (e.g. Kenyon & Luu 1999). 
However, it is expected that growth timescales should increase rapidly with heliocentric 
distance, perhaps as t oc -R^ (Wetherill 1989). One could then expect a reduction in the 
number of large objects beyond 50 AU, as per (2) above, and a correspondingly steeper size 
distribution, as in (1), at larger heliocentric distances. However, with t oc -R^, the timescales 
for growth at -R = 41 AU (inner edge) and -R = 50 AU (outer edge) are only in the ratio 
1.8:1. In addition, we observe no correlation between size and semimajor axis among the 
Classical KBOs. 

To test scenario (1), we took our untruncated best-fit model and varied the size 
distribution index gout for bodies with semimajor axes a > -Rmax, keeping the KBO mass 
across the -Rmax boundary constant. We then found the minimum gout consistent with our 
null detection beyond -Rmax- This mass-conservation model is very sensitive to the chosen 
minimum body radius Tmin, because for any gout > 4, most of the mass is in the smallest 
bodies (Dohnanyi 1969). The minimum size-distribution index required as a function of 
Tmin appears in Table |T2|. If mass is conserved for the observable range of bodies, rmin = 50 
km, the observed edge cannot be explained by a change in the size distribution unless 
g > 10 (3(t), an unphysically large value. For the conservative case of rmin = 6 km (roughly 
the size of cometary nuclei, Jewitt 1997), the observed edge could only be explained by 
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q > 5.6 (3(t) beyond i?max- We know of no population of bodies with a comparably steep 
size distribution. Thus, we conchidc that the observed edge is unlikely to be solely caused 
by a change in the size distribution beyond i^max- 

A similar procedure was followed for possibility (2). Here again, we took our best-fit 
truncated model and extended it to large heliocentric distances. Then, rmax was varied to 
find the largest value that could explain our null detection beyond i?max; keeping the total 
number density of objects with radii r < Tmax constant. Wc found that rmax < 75 km (3(t) 
was required beyond -Rmax to explain the observed edge. This is a factor ~ 5 smaller radius 
and a factor ~ 150 less volume than our largest object found within -Rmax (1999 CDisg, 
~ 400 km in radius). Such a severe change in the maximum object size beyond i?obs would 
have to occur despite the fact that growth timescales vary by less than a factor of ~ 2 over 
the observed Classical KBO range, as explained above. 

One might also expect (3) to be true, as KBO surfaces could darken over time with 
occasional resurfacing by collisions (Luu & Jewitt 1996), and long growth timescales indicate 
long collision timescales as well. However, the geometric albedo would have to be p < 0.008, 
a factor 5 lower than that of the CKBOs in our model, assuming a constant number density 
of objects across the transition region. We are not aware of natural planetary materials 
with such low albedos. 

The dynamical cases (4), a drop in the eccentricity distribution, and (5), a steeper 
ecliptic plane density index, can also be rejected. Even an extreme change in the eccentricity 
distribution cannot explain our observations. Lowering eccentricity from e = 0.15 (a high 
value for the Classical KBOs) to e = results in a perihelion change from 42.5 AU to 50 AU 
for an object with semimajor axis 50 AU. Such a change corresponds to a 0.7 magnitude 
change in perihelion brightness, and to a factor 2.8 change in the surface density of objects 
expected from our a = 0.63 CLF. This model is rejected by our observations at the > 5a 
level. The variation in ecliptic plane surface density with respect to heliocentric distance 
was assumed to follow a power law with index p = 2 in our model. However, even a large 
increase to p = 5 would result in a reduction in surface density of a factor 2.7 in the 41 AU 
to 50 AU range, which can also be rejected as the cause of our observed edge at the > 5a 
level. 

Since scenarios (1) through (5) seem implausible at best, we conclude that the most 
probable explanation for the lack of objects discovered beyond Rmax is (6), the existence of 
a real, physical decrease in object number density. There have been few works considering 
mechanisms for such truncation. The 2:1 mean-motion Neptune resonance at a ~ 47.8 AU 
is quite close to the observed outer edge of the belt. However, given the Neptune resonance 
sweeping model (Malhotra 1995), the resonance could not cause an edge. The sweeping 
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theory predicts that the 2:1 resonance should have passed through the Classical Kuiper 
Belt as Neptune's orbit migrated outwards to its present semimajor axis. Thus, the KBOs 
interior to the current 2:1 resonance (a ~ 47.8 AU) could have been affected by this process, 
but an edge at -Rmax cannot be explained by such a model. Ida, Larwood & Burkert (2000) 
simulate the effect of a close stellar encounter on the Kuiper Belt, suggesting that KBO 
orbits beyond 0.25-0.3 times the stellar perihelion distance would be disrupted and ejected 
for a variety of encounter inclinations. Thus, an encounter with a solar mass star with 
perihelion at ~ 200 AU might explain the observed edge. Such encounters are implausible 
in the present solar environment but might have been more common if the sun formed with 
other stars in a dense cluster. 



6. Constraints on a Distant Primordial Kuiper Belt 

While our observations indicate a dearth of objects beyond 50 AU, it is also possible 
that a "wall" of enhanced number density exists at some large R ^ 100 AU distance, 
as suggested by Stern (1995). We know that the Kuiper Belt has lost much mass since 
formation because the present mass is too small to allow the observed objects to grow in 
the age of the solar system. Kenyon & Luu (1999) found that the primordial Kuiper Belt 
mass in the 30 AU < i? < 50 AU region could have been some ~ 10M@, compared to the 
~ O.IM® we see today. Stern (1995) also speculated that the primordial surface density 
may be present at large heliocentric distances. We model this primordial belt as analogous 
to the CKBOs in terms of eccentricity, inclination and size distribution, but containing a 
factor of 100 more objects and mass per unit volume of space. These objects would be 
readily distinguishable from the rest of the objects in our sample as they would have low 
eccentricities characteristic of the CKBOs (e < 0.25) yet would have very large semimajor 
axes (a > 90 AU). Since we have discovered no such "primordial" objects, Poisson statistics 
state that the 3a upper limit to the sky-area number density of primordial KBOs is 5.9 
in 37.2 sq deg, or 0.16 primordial KBOs dcg~^. We constrain the primordial KBOs by 
allowing the inner edge of the population, amin, to vary outwards, while keeping the outer 
edge fixed at 250 AU. We find that amin = 130 AU coincides with the 3a limit on the inner 
edge of the belt, nearly at the extreme distance limit of our survey. An object discovered 
at our survey magnitude limit mij5o = 23.7 at this distance would have diameter D f» 1800 
km (approximately 25% smaller than Pluto) assuming a 4% albedo. 
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7. Summciry 

New measurements of the Kuiper Belt using the world's largest CCD mosaic array 
provide the following results in the context of our Classical KBO model. 

(1) The slope of the differential size distribution, assumed to be a power law, is 
q = 4.0lo;5 (l*^)- This is consistent with accretion models of the Kuiper Belt (Kenyon 
& Luu 1999). This distribution implies that the surface area, the corresponding optical 
reflected light and thermal emission are dominated by the smallest bodies. 

(2) The Classical KBOs inhabit a thick disk with Half- Width 20°+^o {la). 

(3) The Classical KBOs have a velocity dispersion of l.lGloiie km/s. 

(4) The population of Classical KBOs larger than 100 km in diameter 
^''cKBOs(^ > 100 km) = 3.8tli x 10^ (la). The corresponding total mass of bodies 
with diameters between 100 km and 2000 km is M(100 km, 2000 km) ~ O.O3M0 , assuming 
geometric red albedo pr = 0.04 and bulk density p = 1000 kg m~^. 

(5) The approximate population ratios of the Classical, Scattered, 3:2 Resonant 
(Plutinos) and 2:1 Resonant KBOs are 1.0:0.8:0.04:0.07. 

(6) The Classical Kuiper Belt has an outer edge at i? = 50 AU. This edge is unlikely 
to be due to a change in the physical properties of the CKBOs (albedo, maximum object 
size, or size distribution). The edge is more likely a real, physical depletion in the number 
of bodies beyond ~ 50 AU. 

(7) There is no evidence of a primordial (factor 100 density increase) Kuiper Belt out 
to heliocentric distance R = 130 AU. 
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Fig. 1. — Fields imaged in Feb 1999. The ecliptic is denoted by a solid line. 



-25- 




Fig. 2. — Fields imaged in Sep 1999. The ecliptic is denoted by a solid line. 
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Fig. 3. — Fields imaged in Mar 2000. The ecliptic is denoted by a solid line. 
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Fig. 4. — Inclination vs. discovery distance of all multi-opposition KBOs. The hollow 
circles represent quantities determined from < 90 day tiniebase during the first opposition. 
The connected filled circles represent the orbital solution including second opposition 
observations. Note that for all objects except one, quantities are well determined during 
the first opposition. 
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Fig. 5. — Eccentricity vs. scmimajor axis of all multi-opposition KBOs with a < 50 AU. 
The hollow circles represent the orbits determined during the first opposition. The connected 
filled circles show the orbital elements computed including second opposition observations. 
2 CKBOs were reclassified as Scattered KBOs, and 1 Scattered KBO was reclassified as a 
CKBO. In addition, 3 Resonant KBOs were reclassified as non-resonant objects. 
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Fig. 6. — Eccentricity vs. semimajor axis of all KBOs discovered in this work with semimajor 
axes a < 50 AU. Note that few objects were found in the 3:2 resonance compared to previous 
studies. The area enclosed by a solid line indicates our criteria for selecting Classical KBOs, 
semimajor axes 40.5 AU < a < 46 AU and perihelia > 37 AU. 
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Apporent Red Mognitude 



Fig. 7. — Our measurement of the Cumulative Luminosity Function (CLF), which represents 
the number of KBOs deg^^ near the echptic (filled circles) brighter than a given apparent 
red magnitude. Other points are previous works (see text for abbreviations), with arrows 
denoting upper limits. The line represents a fit to our data alone, yielding a — 0.63 ± 0.06, 
corresponding to g = 4.15 ± 0.3 assuming the the albedo and heliocentric distance 
distributions are independent of the size distribution. Abbreviations are as follows: OOSJTBA 
is Sheppard et al. (2000), 99CB is Chiang & Brown (1999), 98GKNLB is Gladman et al. 
(1998), 98JLT is Jewitt, Luu & Trujillo (1998), 98LJ is Luu & Jewitt (1998), 98TJ is Trujillo 
& Jewitt (1998), 96JLC is Jewitt, Luu & Chen (1996), 95ITZ is Irwin, Tremaine & Zytkow 
(1995), 90LD is Levison & Duncan (1990), 89K is Kowal (1989), 88LJ is Luu & Jewitt 
(1988), and 61T is Tombaugh (1961). 
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Proboblllty Contours 




Slope Porometer q 



Fig. 8. — The maximum likelihood simulation of the size distribution power-law exponent. 
Contours of constant likelihood (la, 2(j, ... 5a) are shown for a model with differential 
size distribution q (x-axis) and total number of objects greater than 100 km in diameter 
N{D > 100 km) (y-axis). The maximum likehhood parameters (denoted by an x) occur at 
q = 4.0 and A^ckbos(^ > 100 km) = 3.8 x lOl 
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Fig. 9. — The Differential Luminosity Function (DLF), equal to the nunAer of KBOs 
deg^^ near the ecliptic (filled circles). Three different models of the observed magnitude 
distribution are plotted from our maximum likelihood model (lines), representing the 
expected DLF for the +la (dotted), best-fit (sohd), and —la (dotted) cases of g = 3.5, 
4.0, and 4.6, respectively. 
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Fig. 10. — Inclination vs. semimajor axis of all KBOs discovered in this work with semimajor 
axes a < 50 AU. 
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1 10 

Object Inclinotion [deg] 

Fig. 11. — The fraction / of an orbit spent within ihl° (solid line), 10° ± 1° (short dashed 
line), and 20° ± 1° (long dashed line) of the ecliptic, as a function of object inclination i. 
The dotted line has a slope of 
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Fig. 12. — The maximum likelihood simulation. Contours of constant likelihood (la, 2ct, 
... 5o") are shown for a model with Gaussian half-width (x-axis) and total number of 
CKBOs with diameters greater than 100 km iVcKBOs(-D > 100 km) (y-axis). The maximum 
hkelihood occurs at iVcKBOs(-D > 100 km) = 3.8 x 10^ and ii/2 = 20°. 
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Fig. 13. — Surface density of KBOs brighter than m^i = 23.7 vs. echptic latitude. The 

+6° 

sohd hne represents the best fit ii/2 = 20°_^o CKBO model while the dotted lines represent 
the la errors. The CKBO model has been multiphed by the observed KBO/CKBO ratio 
(86/49 = 1.76) for display purposes, to simulate the surface density of the more numerous 
KBOs. 
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Fig. 14. — Observed heliocentric discovery distance (data points) and expected discoveries 
assuming the best-fit untruncated CKBO model (solid line). Note the very sharp drop in 
discovery statistics beginning at ~ 46 AU, violating the model. This is consistent with an 
outer edge to the Classical Kuiper Belt at 50 AU (3cr). 
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Table 1. CFHT Survey Parameters 



Quantity 



CFHT 3.6m 



Focal Ratio 
Instrument 
Plate Scale [arc sec/pixel] 
North-South extent [deg] 
East- West extent [deg] 
Field Area [deg^] 
Total Area [deg^] 
Integration Time [sec] 
Readout Time [sec] 

^^[arc sec] 
Filter 
Quantum Efficiency 



f/4 



CFHT 12k X 8k Mosaic 
0.206 
0.47 
0.70 
0.330 

73 
180 

60 
23.7 
0.7-1.1 

R 
0.75 



^the red magnitude at which detection efficiency 
reaches half of the maximum efficiency 

^the typical Full Width at Half Maximum of stellar 
sources for the surveys. 
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Table 2. CFHT 12k Photometric Calibration 
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09 
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8 



Note. — Zero points were consistent 
between observing runs, however, 3 chips were 
replaced and several of the remaining chips 
were shifted in position and renumbered after 
the Feb 1999 run. 

^New chip added after Feb 1999 run. 
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Table 3. CFHT Field Centers 
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Table 3 — Continued 
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Table 3 — Continued 
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09:40 10:37 12:51 


-10 


10:11:50 


01:08:57 


g 






477U290 


iqqq FpK 


1 9 

XZ< 


1)9:45 iU:4z 12:od 


1 n 

-10 


10:15:40 


nn . /i 'v. 1 1 

00:47:11 


g 






477030o 


1 qqq FpK 


1 9 

xz« 


nn en 1 n Ar* 1 nn 

09:50 10:46 13:00 


-10 


10:19:31 


nn c on 

00:25:29 


g 






477(J31o 


iqqq FpK 


1 9 

XZ/ 


nncc in CI 10 nc 

09:55 10:51 13:05 


1 n 

-10 


1 n 00 1 n 

10:23:19 


nn no /i 

00:03:42 


g 






477Uozo 


iqqq Fph 


1 9 

X^ 


nn cn in C£? 10 10 

09:59 10:55 13:lo 


1 n 

-10 


1 n o'? AT 

10:27:07 


nn 1 1 n 

-00:18:10 


g 






477Uooo 


iqqq FpH 


1 9 

X^ 


1 n n /I 11 nn 10 00 

10:04 11:00 13:23 


1 n 

-10 


1 n on c c 

10:30:55 


nn An 01 

-00:40:21 


g 






A ^Tr^o A „ 

477Uo4o 


iqqq Fph 


1 9 

XZi 


1 n nn 1 1 nc 1 o ot 

10:09 11:05 13:27 


1 n 

-10 


1 n /I /1 

10:34:43 


ni no AIT 

-01:02:47 


g 






477U350 


iqqq Fph 


1 9 

XZ< 


inio 11 in io oo 

10:13 11:10 13:32 


1 n 

-10 


in 00 on 

10:38:29 


ni c no 

-01:25:03 


g 






477U3DO 


1 qqq Fph 


1 9 

XZ< 


inio 11 1/1 loo/^ 

10:18 11:14 13:35 


1 n 

-10 


in /I 1 

10:42:15 


n 1 A 'T A 

-01:47:48 


g 






4771590 


iqqq Fph 


1 5 


n/^ CO no 1 n nn nn 

06:53 08:12 09:09 


10 


no in /in 

08:10:49 


on no 01 

30:02:21 


P 






477I6IO 


1999 Feb 


15 


c\i'7 n/^ no 1 f\f\ 1 

07:06 08:17 09:14 


1 n 

10 


no 1 c no 

08:15:02 


on A Cl At 

29:48:41 


P 






4771DZO 


1999 Feb 


15 


U7:ll Uo:zl 09:19 


10 


no, 1 n. 1 

08:19:10 


on.o C , 1 

z9:oo:18 


m 






477IOOO 


1999 Feb 


15 


n'? 1 ^7 no c\£* nn oo 

07:17 08:25 09:23 


1 n 

10 


no 00 OT 

08:23:27 


on on c c 

29:20:55 


P 






4771d4o 


1999 Feb 


15 


n^ o /I no 01 nn oo 

07:24 08:31 09:28 


1 n 

10 


no 0^ 00 

08:27:38 


on n/^ A n 

29:05:45 


P 






4771d5o 


1999 Feb 


15 


f\ f-7 on no 0/^ nn 00 

07:29 08:35 09:33 


1 n 

10 


no 01 /in 

08:31:49 


00 C 1 CI 

28:51:51 


g 






477I660 


1999 Feb 


15 


07:33 08:41 09:37 


10 


08:35:59 


28:36:50 
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4771670 


1999 Feb 


15 


07:38 08:46 09:42 


10 


08:40:07 


28:21:24 


m 






477I680 


1999 Feb 


15 


07:43 08:50 09:46 


10 


08:44:16 


28:06:06 
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4771690 


1999 Feb 


15 


07:50 08:55 09:54 


10 


08:48:23 
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4771 70o 


1999 Feb 


15 


07:56 09:00 09:59 


10 


08:52:29 


27:33:09 


P 
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00 


4771710 


1999 Feb 


15 


08:01 09:04 10:03 


10 


08:56:34 


27:16:38 
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4771980 
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10:17 11:10 12:01 


-10 


09:16:50 


05:48:46 


m 






4771990 


1999 Feb 15 


10:22 11:15 12:06 


-10 


09:20:50 


05:30:28 


g 






4772OO0 


1999 Feb 15 


10:27 11:20 12:11 


-10 


09:24:50 


05:11:47 
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10y<1 1y<OM1CO^ 
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13:34:26 
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1999 I'eb 15 


10 AO t A A^ 1C A'K 
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-04:34:07 
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21:36:32 


-04:15:09 
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10 


21:40:17 
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1999 Sep Ub 
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00:35:33 
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04 
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1999 Sep U6 


09:43 11:00 11:54 
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UU:U9:U7 


00:58:31 
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502100O 


1999 Sep 06 


09:48 11:05 11:59 
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00:12:43 


01:21:56 
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1999 Sep 06 


09:53 10:49 12:03 
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10 
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1999 hep 07 
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1999 Sep 07 


09:54 10:43 11:31 


0.0 


22:35:33 


-08:52:51 


g 
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02 
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1999 Sep 07 


09:59 10:48 11:35 
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22:39:11 


-08:31:00 


m 
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1999 Sep 07 


10:04 10:53 11:40 


0.0 


22:42:48 


-08:09:16 


g 
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1999 Sep 07 


10:09 10:58 11:45 
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22:46:25 
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g 
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10 
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g 
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10 
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01:28:59 


m 


1999 RJ215 


00 


502222O 


1999 Sep 07 


10:28 11:17 12:04 


10 


22:42:29 


01:50:25 


g 
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Table 3 — Continued 



ID 


UT date 


UT times 


/3 


a 




e 


Objects 


Cliip 


5022260 


1999 Sep 07 


10:33 11:21 12:09 


10 


oo An o/^ 

22-Ab.2b 


AO 1 O O/^ 

02:12:36 


g 






o(Jzz44o 


1999 


Sep 


07 


12:15 


13:11 


14:10 


-1(J 


AA. yl C 1 '7 

(JU:45:17 


A C . AO . OO 


m 


1 A A A T~> T7" 

1999 KK215 


AO 

02 
























1 f\r\f\ Ti AT 

1999 RN215 


11 


rrnOO /I FC^ 

o(Jzz4oo 


1999 


Sep 


07 


12:20 


13:26 


14:15 


-iU 


AA. /I Q. 


-U4:4o:Uo 


m 






5Uz24do 


1999 


Sep 


07 


12:25 


13:31 


14:20 


-10 


A A CO OC 

00:52:35 


A A OO OA 

-04:22:30 


P 






5Uz247o 


1999 


Sep 


07 


12:29 


13:36 


14:25 


-10 


00:56:15 


-03:59:34 


m 






5U22480 


1999 


Sep 


07 


12:34 


13:41 


14:30 


0.0 


AA A r" 1 

00:45:16 


r\ A r 1 A r\ 

04:51:49 


m 


1 AAA T> "r> 

1999 KK215 


A A 

09 


5U2z49o 


1999 


Sep 


07 


12:39 


13:46 


14:34 


0.0 


AA AO r C 

00:48:55 


AT 1/1 r o 

05:14:52 


m 


1 AAA T^T^ 

1999 KI215 


AO 

02 


b()22b()o 


1999 


Sep 


07 


12:44 


13:51 


14:39 


0.0 


A A . C A . O /I 

00:52:34 


A C . OT. C A 

05:37:50 


m 


1 AAA T~) T T 

1999 KU215 


04 




1999 


Sep 


07 


12:48 


13:55 


14:44 


10 


A A . /I C 1 

00:45:16 


1 /I . C 1 . /I A 

14:51:49 


g 


1 AAA T~>\ T 

1999 KV215 


A 1 

01 
























1999 KW215 


no 
03 


oUziozo 


1999 


Sep 


07 


12:53 


14:00 


14:49 


10 


AA AO C C 

00:48:55 


1 C 1 /I o^ 

15:14:26 


m 






oUzzooo 


1999 


Sep 


07 


12:59 


14:05 


14:54 


10 


A A . C A . O C 

00:52:35 


1 C . 0*7. C A 

15:37:50 


m 






5023670 


1999 


Sep 


08 


06:19 


07:45 


08:55 


-10 


21:25:14 


or" r~\r\ rr^ 

-25:09:59 


P 






502368o 


1999 


Sep 


08 


06:23 


07:50 


09:00 


-10 


o 1 c\f\ r\r\ 

21:29:00 


O /I f 'I O /I 

-24:51:34 


P 






5023710 


1999 


Sep 


08 


06:43 


07:55 


09:05 


-10 


01 OO Ai"? 

21:32:47 


C% A O O O^ 

-24:33:37 


P 






OUZO i zo 


1999 


Sep 


08 


06:48 


08:00 


09:10 


-lU 


z1:od:oz 


-z4:lo:lo 


P 






502374O 


1999 


Sep 


08 


06:59 


08:15 


09:30 


-0.5 


21:36:32 


1 A . AK .f^ A 

-14:45:24 


m 






502375O 


1999 


Sep 


08 


07:04 


08:20 


09:35 


-0.5 


21:40:17 


1 A O/^ O^ 

-14:26:37 


P 


1 AAA T>"V" 

1999 KA215 


A^ 

07 


50237DO 


1999 


Sep 


08 


07:09 


08:25 


09:40 


A CT 

-0.5 


1 A'T A^ 

21:47:45 


1 O /I O O /I 

-13:48:24 


P 






0023/ 


1999 


Sep 


08 


07:14 


08:30 


09:46 


-0.5 


2i:ol:29 


1 Q .on. 1 /I 
-13:29:14 


P 


1999 KY215 


01 
























1999 KZ253 


no 
03 


50237OO 


1999 


Sep 


08 


07:19 


08:34 


09:54 


-10 


21:40:17 


-23:56:27 


P 






5023790 


1999 


Sep 


08 


07:23 


08:39 


09:59 


-10 


21:44:01 


00.0'7.00 

-26.6i:62 


P 






5023o0o 


1999 


Sep 


08 


07:28 


08:43 


10:03 


-10 


01 . AT. AK 

21:47:45 


OO .1 O . O /I 

-zd:lo:z4 


P 


1999 KZ215 


06 


c! no o o 1 ^ 

oUzoolo 


1999 


Sep 


08 


07:33 


08:51 


10:08 


1 A 

-10 


21:51:28 


OO . C A. A /I 

-22:59:04 


P 






502422o 


1999 


Sep 


08 


11:21 


12:23 


13:30 


9.5 


23:33:18 


06:36:36 


P 






502423O 


1999 


Sep 


08 


11:26 


12:28 


13:35 


9.5 


23:36:53 


06:59:39 


m 






502424O 


1999 


Sep 


08 


11:31 


12:32 


13:40 


9.5 


23:40:27 


07:23:07 


P 






502425O 


1999 


Sep 


08 


11:36 


12:37 


13:44 


0.5 


23:33:13 


-02:23:38 


P 


1999 RA216 


09 


502426O 


1999 


Sep 


08 


11:41 


12:42 


13:47 


0.5 


23:36:48 


-02:00:38 


m 






502427O 


1999 


Sep 


08 


11:45 


12:46 


13:52 


0.5 


23:40:24 


-01:37:25 


P 


1999 RB216 


07 
























1999 RC216 


11 


502428O 


1999 


Sep 


08 


11:53 


12:51 


13:57 


9.5 


23:44:04 


07:46:06 


P 






502429O 


1999 


Sep 


08 


11:58 


13:08 


14:50 


9.5 


23:47:39 


08:09:21 


P 






502430O 


1999 


Sep 


08 


12:04 


13:13 


14:55 


9.5 


23:51:13 


08:32:54 


P 
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Table 3 — Continued 



ID 


UT date 


UT times 


(3 ^ 


a ° 


S ^ 


9 °- 


Objects 


Chip 




1QQQ Son 


1 9-nQ 1^-17 1 ii-nn 

-LZ.Ui? -LO.±l 10. UU 


y.u 


zo.O'^.^y 


UO.UU. 14 


P 








1QQQ Qpn 


19-1^ 1 ^-09 1 'i-n'i 
j.z.±o 10. zz ±o.uo 


y.o 


zo.oo.zo 


nO-l Q-94 
uy . ly .z^ 


P 








1QQQ CIpn 
iyyy ot;p uo 


19-18 1 ^-9(1 1 ii-1 n 

±Z.J.O -LO.ZU -LU.-LU 


y.u 


nn. ni .KQ 
uu.ui.uy 


uy.4:Z.o / 


P 






OZ ( i ( 00 


zuuu ivicir oi 


nQ-/i9 1 n-/io 1 1 -"^o 
uy.4z 1U.4Z li. OZ 


on 
zu 


1 1 .CI .OR 
ll.Ol.UD 


on-^7-i^^ 
zu.o ( .00 


g 






11971 74n 


^UUU iVlOji Oi- 


nQ-^'^ 1 n-4fi 1 1 •'^7 

Uy.tJO i-U.^U X-L.O/ 


n n 

u.u 


11-^1 -nfi 

±±.o±.uu 


UU.(J t .04: 


or 
6 


^uuu r -A.53 


05 
uo 


^971 7^n 
(JZ tit ou 


zuuu IVldl 1 


Uc/.oy lu.oz iL.^i 


n n 

u.u 


11. uo.ou 


nn-zLo-1 

UU.4Z. lU 


in 






111 UL) 


9nnn Mar ^1 

zuuu iVldil O-L 


1 n-n4 1 n- lift 11 -a^ 

-LU.U^ -LU.UU ±J..^U 


n n 

u.u 


1 1 •'\^-^A 


UU.ZU.OO 


P 






c;97'?ni r> 
OZ / ou±o 


9nnn At->t no 
zuuu Apr uz 


ns-/ic; nQ-/in ^^^■At^ 
uo.^to uy.'iu lu.^io 


on 
zu 


1 o-nn-/i^ 

IZ.UU.^to 


ly. 00. OZ 


m 






OZ ( ouzo 


onnn At-.t no 
zuuu iipr uz 


ns-c;n nQ-/ifi in-f;n 
uo.ou uy.^io iu.ou 


on 
zu 


1 o-n^-n7 

IZ.UO.U 1 


iy.oy.4o 


m 






OZ / oUoO 


onnn Anv n9 
zuuu Apr uz 


ns-i^fi nQ-f;i in-^i^ 

Uo.OD uy.oi iU.OO 


on 
zu 


IZ.UO.Ol 


1 Q- 0/1 - no 
ly .z^t.uy 


m 






tJZ / OU4:U 


9nnn Anr n9 

zuuu Api uz 


nQ-n9 no-i^^ 1 1 -nn 

uy.uz uy.oo ii.uu 


on 
zu 


IZ.U 1 .(JO 


ly .uo.oz 


m 






OZ 1 oUOO 


onnn A t-vt- no 
zuuu Apr uz 


no-nsj in-ni 
uy.Uo iu.ui ii.uo 


n n 
u.u 


1 o-nn-/i Q 

1Z.UU.40 


nn-n/i -^7 

-UU.U4.0 / 


g 






OZ ( oUDO 


onnn Am- no 
zuuu Apr uz 


nQ-1 "i 1 n-nR 1 1 -i n 
uy.io lu.uo li.iu 


n n 
u.u 


1 o-n'^-n7 

IZ.UO.U 1 


nn-on-1 c; 
-uu.zu.io 


m 


iyy4 vjVg 


1 1 
11 


OZ / oU ( 


onnn At->t no 
zuuu Apr uz 


no-i in-iQ 11-ic; 
uy.iy iu.iy 11.10 


n n 
u.u 


1 o-nc;-'?n 

IZ. UO.OU 


-uu.00.01 


in 


onnn n'w. .» 

zuuu Vjtx\.i47 


1 1 
11 


OZ / oUoO 


onnn At->v no 
zuuu Apr uz 


nQ-9/i 1 n-os 1 1 -01 
uy.z^i lu.zo 11. zi 


n n 
u.u 


1 o-n7-c;i^ 
iz.u / .00 


nn-c;i -os 
-uu.oi.zo 


P 








9nnn Anr 02 


0Q-S1 in-3S 1 1 -97 

Uty.tji xvj.tjo J.x.^( 


90 


12:10:18 


18-52-57 


Ti 
P 






5273100 


2000 Apr 02 


09:36 10:38 11:33 


20 


12:12:43 


18:37:21 


P 






5274550 


2000 Apr 03 


08:31 09:21 10:12 


20 


12:31:55 


16:33:22 


m 


2000 GM147 


03 


5274560 


2000 Apr 03 


08:36 09:26 10:17 


20 


12:34:19 


16:17:59 


g 






5274570 


2000 Apr 03 


08:41 09:31 10:23 


20 


12:36:43 


16:02:37 


g 






5274580 


2000 Apr 03 


08:46 09:36 10:34 


0.0 


12:31:55 


-03:26:38 


g 


2000 GWi46 


09 


5274590 


2000 Apr 03 


08:51 09:41 10:39 


0.0 


12:34:19 


-03:42:01 


m 


2000 GYi46 


04 
















2000 GXi46 


06 


527460O 


2000 Apr 03 


08:56 09:46 10:44 


0.0 


12:36:43 


-03:57:23 


m 






5274610 


2000 Apr 03 


09:01 09:51 10:49 


0.0 


12:39:07 


-04:12:43 


g 






5274620 


2000 Apr 03 


09:06 09:57 10:55 


20 


12:39:07 


15:47:17 


P 






5274630 


2000 Apr 03 


09:11 10:02 11:00 


20 


12:41:31 


15:31:59 


g 






5274640 


2000 Apr 03 


09:16 10:07 11:05 


20 


12:43:56 


15:16:43 


g 
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Table 3 — Continued 



ID 


UT date 


UT times 


l3 


a ^ 






Objects Chip 


5274870 


2000 Apr 03 


11:22 


12:16 13:16 


0.0 


13:54:20 


-11:43:09 


m 




5274880 


2000 Apr 03 


11:26 


12:22 13:21 


0.0 


13:56:48 


-11:56:34 


m 




5274890 


2000 Apr 03 


11:32 


12:27 13:26 


20 


13:54:20 


08:16:52 


m 




527500O 


2000 Apr 03 


12:32 


13:32 14:10 


20 


13:56:48 


08:03:26 


m 




5274910 


2000 Apr 03 


11:45 


12:37 13:37 


20 


13:59:16 


07:50:04 


g 




5274920 


2000 Apr 03 


11:50 


12:42 13:42 


20 


14:01:44 


07:36:48 


m 




5274930 


2000 Apr 03 


11:55 


12:48 13:47 


20 


14:04:13 


07:23:38 


m 




5274940 


2000 Apr 03 


12:00 


12:53 13:52 


20 


14:06:41 


07:10:32 


m 




5274950 


2000 Apr 03 


12:05 


12:59 13:58 


0.0 


14:04:13 


-12:36:23 


P 


2000 GZi46 11 


5274960 


2000 Apr 03 


12:11 


13:04 14:03 


0.0 


14:06:41 


-12:49:28 


P 





Note. — This table lists fields imaged with the CFHT 12k Mosaic camera. Fields were imaged 
in triplets, with UT times given for each image. KBOs found are listed after the field of discovery. 
If more than one KBO was found in each field, they are listed on successive lines. The complete 
version of this table is in the electronic edition of the Journal. The printed edition contains only a 
sample. 

'^J2000 ecliptic latitude, degrees 
'^J2000 right ascension, hours 
"^J2000 declination, degrees 

*^Seeing category: g, m, p represent the good (< 0.8 arc sec), medium (> 0.8 arc sec and < 1.0 arc 
sec), and poor (> 1.0 arc sec) seeing cases, respectively. The efficiency functions for each of these 
cases are presented in Table ^ 
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Table 4. CFHT Survey Efficiency 



Quantity 


Good 


Medium 


Poor 


Global 


Median PSF FWHM ["] 


0.76 


0.90 


1.07 


0.84 


PSF FWHM Range ["] 


0.56-0.80 


0.80-0.99 


1.00-1.40 


0.56-1.40 




0.83 


0.83 


0.83 


0.83 




24.01 


23.64 


23.35 


23.74 


a 


0.29 


0.38 


0.47 


0.48 


Fields Imaged 


95 


89 


49 


233 
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Table 5. CFHT Discovery Conditions 



code 




Date 




R 


A 


a' 


ruR ± a 


mR(l,l,0) 


D 


MPC name 




[deg] 








[AU] 


[AU] 


[deg] 






[km] 




C17000 


+10 


1999 


02 


15 


41.161 


42.085 


0.5 


21.36 + 0.15 


5.16 


527 


1999 CCi58 


C71710 


0.0 


1999 


02 


10 


47.562 


48.474 


0.4 


21.08 + 0.17 


4.27 


797 


1999 CDi58 


C72105 


0.0 


1999 


02 


10 


37.694 


38.631 


0.5 


22.80 ± 0.09 


6.99 


228 


1999 CViis 


C72505 


0.0 


1999 


02 


10 


46.521 


47.477 


0.3 


23.54 ± 0.32 


6.82 


246 


1999 CPi53 


C72600 


0.0 


1999 


02 


10 


37.713 


38.671 


0.4 


23.76 ± 0.14 


7.94 


146 


1999 CKi58 


C75803 


0.0 


1999 


02 


10 


34.141 


35.053 


0.6 


23.61 + 0.13 


8.22 


129 


1999 CY118 


C75905 


0.0 


1999 


02 


10 


45.106 


46.010 


0.5 


24.40 ± 0.27 


7.82 


155 


1999 CZ118 


C76103 


0.0 


1999 


02 


10 


44.262 


45.153 


0.5 


24.11 + 0.26 


7.61 


171 


1999 CA119 


C76104 


0.0 


1999 


02 


10 


42.565 


43.455 


0.6 


24.05 + 0.21 


7.72 


162 


1999 CW118 


C76107 


0.0 


1999 


02 


10 


40.298 


41.190 


0.6 


22.65 ± 0.07 


6.55 


278 


1999 CB119 


C76207 


0.0 


1999 


02 


10 


43.531 


44.415 


0.6 


22.97 ± 0.06 


6.54 


280 


1999 CC119 


C76502 


0.0 


1999 


02 


10 


44.301 


45.161 


0.6 


23.56 ± 0.14 


7.06 


220 


1999 CD119 


C76609 


0.0 


1999 


02 


10 


44.002 


44.852 


0.6 


24.11 + 0.11 


7.64 


169 


1999 CX118 


C76800 


0.0 


1999 


02 


10 


28.873 


29.711 


1.0 


23.19 + 0.04 


8.52 


112 


1999 CE119 


C76907 


0.0 


1999 


02 


10 


42.953 


43.778 


0.7 


23.96 ± 0.08 


7.59 


172 


1999 CW131 


C79710 


+0.5 


1999 


02 


10 


— 


— 


— 


23.35 + 0.14 


— 


— 


lost 


C79900 


+0.5 


1999 


02 


10 


— 


— 


— 


22.78 ± 0.12 


— 


— 


lost 


C85003 


+0.5 


1999 


02 


11 


— 


— 


— 


23.68 + 0.13 


— 


— 


lost 


C85110 


+0.5 


1999 


02 


10 


38.774 


39.705 


0.5 


22.75 + 0.20 


6.82 


246 


1999 CF119 


C85200 


+0.5 


1999 


02 


10 


41.806 


42.739 


0.4 


23.63 + 0.07 


7.37 


191 


1999 CG119 


C85204 


+0.5 


1999 


02 


10 


32.938 


33.875 


0.5 


21.85 + 0.09 


6.61 


271 


1999 CLi58 


C85300 


+0.5 


1999 


02 


10 


45.906 


46.845 


0.4 


23.88 ± 0.23 


7.22 


204 


1999 CH119 


C85404 


+0.5 


1999 


02 


10 


42.175 


43.122 


0.4 


23.98 ± 0.08 


7.68 


165 


1999 CS153 


C85410 


+0.5 


1999 


02 


11 


— 


— 


— 


23.82 ± 0.08 


— 


— 


lost 


C85504 


+0.5 


1999 


02 


10 


41.130 


42.082 


0.4 


24.24 + 0.16 


8.05 


139 


1999 CK119 


C85509 


+0.5 


1999 


02 


10 


41.450 


42.400 


0.4 


23.03 + 0.12 


6.80 


247 


1999 CJ119 


C85600 


+0.5 


1999 


02 


10 


45.773 


46.727 


0.3 


22.52 ± 0.03 


5.87 


381 


1999 CLu9 


C85700 


+0.5 


1999 


02 


10 


41.159 


42.117 


0.3 


23.25 + 0.10 


7.05 


220 


1999 CM119 


C85808 


+0.5 


1999 


02 


10 


45.122 


46.085 


0.3 


23.41 ± 0.12 


6.82 


246 


1999 CG154 


C85904 


+0.5 


1999 


02 


10 


43.697 


44.664 


0.3 


23.96 + 0.24 


7.51 


179 


1999 CN119 


C85907 


+0.5 


1999 


02 


10 


42.168 


43.135 


0.3 


22.96 + 0.10 


6.66 


264 


1999 CX131 


C85909 


+0.5 


1999 


02 


11 








23.99 ± 0.05 






lost 


C88505 


0.0 


1999 


02 


10 


40.502 


41.482 


0.2 


23.11 + 0.21 


6.98 


228 


1999 CM153 


C88600 


0.0 


1999 


02 


10 


37.016 


37.995 


0.2 


23.61 ± 0.07 


7.87 


151 


1999 CYiai 


C88608 


0.0 


1999 


02 


11 








23.64 ± 0.33 






lost 


C88902 


0.0 


1999 


02 


10 


41.998 


42.968 


0.2 


23.99 ± 0.08 


7.71 


163 


1999 CZiai 
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Table 5 — Continued 



code 


/5 
[dcg] 


Date 




R 
[AU] 


A 

[AU] 


a' 
[deg] 


ruR ± a 




mR(l,l,0) 


D 
[km] 


MPC 


I name 


C88905 


0.0 


1 QQQ 


09 
uz 


11 








24.20 ± 0. 


,28 






lost 




dr\r\r\r\ 

C89000 


0.0 


1 QQQ 




10 


42.306 


43.273 


0.3 


23.75 ± 0. 


1 A 

,14 


7.44 


185 


1999 


CNi53 


C89301 


0.0 


1 QQQ 


09 
uz 


10 


38.704 


39.660 


0.4 


22.91 ± 0. 


,26 


6.97 


229 


1999 


CAi32 


C89500 


0.0 


1 QQQ 




10 


43.383 


A A oor\ 

44.330 


0.4 


23.99 ± 0. 


O A 

,04 


7.57 


1 TO 

173 


1 OOO 

1999 




C89503^ 


0.0 


1 QQQ 




10 


y< C\ At 

44.241 


AT" 1 OT 

45.187 


C\ A 

0.4 


OO C% A 1 O 

23.24 ± 0. 


,11 


^O 

6.73 


256 


1 OO r" 

1995 


DCo2 


C89507 


0.0 


1 QQQ 


09 
uz 


10 


30.805 


31.752 


0.5 


22.11 ± 0. 


,09 


7.15 


211 


1999 




C92411 


-0.5 


1 QQQ 


09 
uz 


10 


21. I'd" 


28.719 


0.7 


o 1 on 1 r\ 

21.89 ± 0. 


,22 


7.38 


190 


1999 


CMi58 


f~\(\Ci r\c\ 

C98602 


-0.5 


1 QQQ 


09 
uz 


12 


43.113 


44.031 


0.5 


oo 1 r\ 

22.96 ± 0. 


,10 


6.57 


276 


1999 




C98603 


-0.5 


IQQQ 


09 
uz 


12 


39.970 


40.889 


0.5 


OO or' 1 o 

22.95 ± 0. 


,17 


OO 

6.88 


ooo 

239 


1 OOO 

1999 


COi53 


C98703 


-0.5 


1 QQQ 


09 


12 


41.324 


42.250 


0.5 


23.85 ± 0. 


,14 


7.64 


168 


1999 




C99411 


-0.5 


1 QQQ 


09 
uz 


12 


41.457 


42.418 


0.3 


oo -t i~\ 1 o 

23.19 ± 0. 


,18 


6.97 


229 


1999 


CUi53 


C99502 


-0.5 


1 QQQ 


09 

uz 


12 


42.006 


AC\ (\ ^(\ 

42.970 


0.3 


OO ^o 1 r\ 

23.73 ± 0. 


,25 


7.45 


184 


1999 


UMi54 


CU47U0 





1 QQQ 


OQ 


Ob 


42.054 


A o c\n A 

42.964 


0.6 


OO O O 1 M 

23.83 ± 0. 


,25 


7.54 


176 


1999 


nD2i4 


C04707 





1 QQQ 


OQ 


06 


39.838 


40.754 


0.6 


o o /I o 1 r\ 

23.49 ± 0. 


1 A 

,14 


7.45 


184 


1999 


i<'-L214 


C04910 





1 QQQ 


OQ 


06 


44.240 


Aft f*^ 

45.167 


0.5 


OO O/^ 1 r\ 

22.96 ± 0. 


OO 

,23 


6.45 


001 

291 


-1 OOO 

1999 


RU214 


C05105 





1 QQQ 


OQ 


06 


47.685 


48.626 


r\ A 

0.4 


O A OO 1 O 

24.08 ± 0. 


OO 

,20 


l'7 O/^ 

7.26 


001 

201 


1 OOO 

1999 


RV214 


C09804 





IQQQ 


OQ 


06 


41.990 


AC\ f\ A A 

42.944 


r\ A 

0.4 


OO Ai'7 1 O 

23.47 ± 0. 


O A 

,34 


^ 1 O 

7.19 


OOT 

207 


1 OOO 

1999 


RW214 


/~\ "i f\t~i r\f\ 

C10300 





1 QQQ 


OQ 


06 


44.694 


45.621 


0.5 


OO OO 1 o 

22.90 ± 0. 


,20 


6.36 


304 


1999 


KA214 


f i -\ r\o r\'\ 

C10801 


+10 


1 QQQ 


OQ 


06 


36.884 




0.7 


22.76 ± 0. 


,08 


7.03 


223 


1999 


RY214 


UlL)9L)8 


+10 


1 QQQ 


OQ 


06 


39.016 


39.906 


0.7 


OO /I o 1 r\ 

26A2 ± 0. 


1 O 

,10 


'T An 


1 O O 

183 


1 OOO 

1999 


KZi2i4 


C13601 





1 QQQ 


OQ 


06 


30.686 


31.610 


0.7 


OO 1 O 

23.71 ± 0. 


OO 

,22 


O ^O 

8.78 


OO 

99 


-1 OOO 

1999 


±<'i>215 


C13604 





1 QQQ 


OQ 


06 


42.465 


43.380 


0.5 


OO OO 1 o 

22.88 ± 0. 


,15 


6.55 


279 


1999 


l^*-^215 


C14209 


+10 


1 QQQ 


OQ 


06 


37.824 


38.692 


0.8 


O O O 1 o 

23.08 ± 0. 


,11 


7.24 


202 


1999 


RD215 


C18dU4 


U 


1 QQQ 


OQ 


07 


41.618 


/I o r o o 

42.588 


n A 

0.4 


OO r '7 1 O 

22.57 ± 0. 


,17 


n O o 

6.32 


309 


1 OOO 

1999 


ID 


C18606 





1999 


09 


07 


/I o r\r\r\ 

43.000 


43.964 


0.4 


o o o r' 1 r\ 

22.95 ± 0. 


,07 


6.57 


276 


1999 


^V'i^215 


C21502 





1999 


09 


07 


44.015 


45.017 


0.2 


O O /I O 1 o 

23.49 ± 0. 


,13 


6.98 


228 


1999 


TD 

i»''j215 


C21903 


+10 


1999 


09 


07 


36.161 


37.154 


0.3 


i%0 O -1 1 

23.91 ± 0. 


,21 


8.27 


126 


1999 


RH215 


C22100 


+10 


1999 


09 


07 


34.047 


35.041 


0.3 


22.42 ± 0. 


,08 


7.03 


223 


1999 


RJ215 


C24402 


-10 


1999 


09 


07 


42.041 


42.952 


0.6 


23.33 ± 0. 


,11 


7.04 


222 


1999 


RK215 


C24411 


-10 


1999 


09 


07 


43.966 


44.869 


0.6 


23.02 ± 0. 


,11 


6.56 


277 


1999 


RN215 


C24809 





1999 


09 


07 


39.498 


40.391 


0.7 


23.93 ± 0. 


,07 


7.91 


148 


1999 


RR215 


C24902 





1999 


09 


07 


42.205 


43.089 


0.6 


23.12 + 0. 


,19 


6.83 


245 


1999 


RT215 


C25004 





1999 


09 


07 


39.515 


40.383 


0.7 


22.87 + 0. 


,20 


6.82 


245 


1999 


RU215 


C25101 


+10 


1999 


09 


07 


35.622 


36.466 


0.9 


23.73 + 0. 


,19 


8.16 


132 


1999 


RV215 


C25103 


+10 


1999 


09 


07 


32.291 


33.139 


1.0 


23.31 + 0. 


,16 


8.17 


132 


1999 


RW215 


C37507 


-0.5 


1999 


09 


08 


41.130 


42.059 


0.5 


23.42 + 0. 


,16 


7.23 


203 


1999 


RX215 
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Table 5 — Continued 



code 


/3 


Date 


R 


A 


a' 


rriR + a 


mfi(l,l,0) 


D 


MPC name 




[deg] 






[AU] 


[AU] 


[deg] 






[km] 




C37701 


-0.5 


1999 


09 08 


36.314 


37.262 


0.5 


22.25 ±0.12 


6.60 


272 


1999 RY215 


C37703 


-0.5 


1999 


09 08 


39.955 


40.905 


0.5 


22.04 ±0.12 


5.97 


364 


1999 RZ253 


C38006 


-10 


1999 


09 08 


30.068 


30.976 


0.8 


22.24 ±0.16 


7.40 


188 


1999 RZ215 


C41409 


+0.5 


1999 


09 08 


42.379 


43.377 


0.2 


23.08 ±0.17 


6.74 


255 


1999 RA216 


C41607 


+0.5 


1999 


09 08 


33.725 


34.718 


0.3 


22.42 ±0.11 


7.05 


221 


1999 RB216 


C41611 


+0.5 


1999 


09 08 


46.992 


47.985 


0.2 


23.56 ±0.11 


6.79 


249 


1999 RC216 


D17405 


0.0 


2000 


03 31 


37.749 


38.722 


0.3 


23.33 ±0.16 


7.51 


179 


2000 FX53 


D30611'' 


0.0 


2000 


04 02 


41.298 


42.278 


0.3 


23.39 ±0.10 


7.18 


208 


1994 GVg 


D30711 


0.0 


2000 


04 02 


33.551 


34.532 


0.3 


23.56 ±0.12 


8.24 


127 


2000 GK147 


D45503 


+20 


2000 


04 03 


37.051 


37.976 


0.6 


22.52 ±0.20 


6.78 


250 


2000 GM147 


D45809 


0.0 


2000 


04 03 


40.209 


41.205 


0.1 


23.81 ±0.15 


7.71 


163 


2000 GWi46 


D45904 


0.0 


2000 


04 03 


43.753 


44.750 


0.1 


23.85 ±0.11 


7.40 


188 


2000 GYi46 


D45906 


0.0 


2000 


04 03 


43.605 


44.602 


0.1 


23.18 ±0.11 


6.74 


255 


2000 GXi46 


D49511 


0.0 


2000 


04 03 


42.935 


43.878 


0.4 


23.36 ±0.15 


6.99 


228 


2000 GZi46 



'^These are previously known objects serendipitously imaged in survey fields. 

Note. — mpt is the red magnitude of the object, with 1 sigma error a. m/j(l,l,0) is the limiting red 
magnitude at geocentric distance A = 1 AU, heliocentric distance R = 1 AU, and phase angle a' = 0, 
computed from discovery geometry. Diameter D is computed directly from m/j( 1,1,0) via Equation ^, 
assuming p = 0.04. Some quantities were not computed for lost objects because observations span only 2 
hours. 
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Table 6. CFHT Orbital Elements 



code 


a 


e 


i 


n 


id 


M 


MJD 




MPC 


! name 


sim*^ 




[AU] 




[dcg] 


[deg] 


[deg] 


[dcg] 












C17000 


54.81897 


0.29808 


18.754 


336.974 


98.551 


29.501 


51800 


2 


1999 






C71710 


44.00641 


0.15285 


25.441 


119.004 


137.939 


236.881 


51800 


2 


1999 


CDi58 


i,q 


C72105 


52.82586 


0.28892 


5.467 


305.621 


1 AC) O 

148.820 


16.449 


51600 


2 


1999 


CViis 




C72505 


45.15557 


0.16213 


3.038 


122.823 


C% Ai^ 1 1^/^ 

247.176 


100.913 


51600 


2 


1999 


CPl53 


i,q 


C72600 


Af\ r\ 1 1 orv 

40.91130 


0.06725 


17.068 


127.720 


322.119 


OP" y< O O 

35.433 


51800 


2 


1999 


CKi58 


i,q 


C75803 


89.81583 


0.61493 


25.619 


163.125 


15.456 


357.553 


51600 


2 


1999 


CYiis 




C75905 


110.36855 


0.65629 


. ■— • C)C\ A 

27.804 


O /I P* 1 '^'7 

345.177 


235.542 


350.265 


51600 


2 


1999 






C76103 


45.15314 


0.0 


0.283 


24.858 


O 1 n o o o 

319.882 


181.758 


51260 


(86) 


1999 




i,q 


C76104 


43.45475 


0.0 


0.819 


1 P" A i^t 

154.471 


1 O 1 

173.481 


1 rvo Trip" 

198.705 


51260 


(86) 


1 f\(\f\ 

1999 


CWiis 


i,q 


C76107 


43.64467 


0.05625 


9.160 


168.036 


358.276 


0.071 


51240 


(62) 


1 f\f\f\ 

1999 


^i>119 


i,q 


C76207 


A A A Adl f 

44.44015 


0.00127 


C\ AV Ci 

0.458 


1 r\r\ o/^ 1 

190.261 


272.993 


65.792 


51800 


2 


1999 




i,q 


C76502 


43.99328 


0.02655 


2.396 


166.381 


183.867 


-ton cx^r" 

180.075 


51240 


(62) 


1999 


CDiig 


i,q 


C76609 


43.66175 


0.02725 


1.784 


175.714 


1 ^P' ^ /i O 

175.748 


180.064 


51240 


(62) 


1999 




i,q 


''vr* a r\f\ 

C76800 


39.33216 


0.24461 


1.429 


171.557 


1.040 


0.083 


51240 


(58) 


1999 


L/Hjlig 




C76907 


/( O -\ t\ A A A 

43.19444 


0.01577 


7.952 


174.627 


209.869 


150.003 


51600 


2 


1 f\f\f\ 

1999 


CWisi 


i,q 


C85110 


91.70394 


0.58203 


19.700 


303.437 


203.636 


o r y< r\i p" 

354.915 


51600 


2 


1999 


119 




C85200 


51.35701 


0.33914 


16.656 


304.264 


257.457 


317.562 


51600 


2 


1 f\f\f\ 

1999 


*-^*J119 




C85204 


41.77052 


0.21889 


10.025 


120.074 


325.711 


26.075 


51800 


2 


1999 


/-IT 

(-^^^158 




C853UU 


43.39088 


0.08942 


19.968 


122.349 


156.415 


O 1 n '1 A'T 

210.747 


CT 1 r\r\ 

51600 


2 


1999 




i,q 


C85404 


44.51143 


0.11926 


0.983 


343.336 


o rj o o o o 

223.233 


293.035 


51600 


2 


1999 




i,q 


C85504 


A o r\o 1 r' p" 

42.08155 


0.0 


11.592 


123.077 


2.999 


0.000 


51220 


(8) 


1 f\f\f\ 

1999 


CKiig 


i,q 


085509 


45.42766 


0.06934 


3.200 


O 1 O O PTV 

313.350 


189.726 


^ Af* C\Ci T" 

346.285 


51600 


2 


1999 


*^Jll9 


i,q 


C85600 


46.95343 


0.02034 


26.292 


125.157 


284.096 


76.274 


51600 


2 


1999 


(-^^^119 




C85700 


44.41239 


0.13382 


2.744 


110 A f\C\ 

118.490 


i~\r\ A /I p" o 

294.453 


61.472 


51600 


2 


1999 


CMiig 


i,q 


C85808 


43.04996 


r\ /~\ o c\c\ 

0.08599 


0.766 


1 r\r\ i^rw 

100.205 


176.323 


o 1 r» r\r\ f 

219.095 


51600 


2 


1999 


ULii54 


i,q 


C85904 


43.89874 


r\ f\'\ '1 A A 

0.01744 


0.997 


347.227 


323.092 


180.915 


51600 


2 


1999 


CNiig 


i,q 


0859U7 


41.48466 


0.18854 


9.760 


127.984 


11/1 o^jn 
114.269 


o^jn TOO 
269.788 


51600 


2 


1999 


/' I V 

*-^^131 




C88505 


44.12661 


0.05994 


0.190 


87.878 


60.139 


0.067 


51240 


(57) 


1999 


CMi53 


i,q 


C88600 


43.92279 


0.13496 


25.167 


148.245 


0.110 


359.999 


51220 


(8) 


1999 


CYisi 


i,q 


C88902 


44.41005 


0.03248 


2.371 


151.799 


359.878 


359.999 


51220 


(8) 


1999 


CZi3i 


i,q 


C89000 


42.59212 


0.01599 


7.480 


153.050 


179.417 


180.056 


51240 


(56) 


1999 


CNi53 


i,q 


C89301 


43.96574 


0.09794 


12.071 


154.871 


0.530 


0.067 


51240 


(56) 


1999 


CAi32 


i,q 


C89500 


44.33005 


0.0 


0.208 


181.177 


336.136 


0.000 


51220 


(30) 


1999 


CQl53 


i,q 


089503^ 


43.86806 


0.06335 


2.349 


154.182 


124.958 


244.781 


51200 


5 


1995 


DCo2 


i,q 


C89507 


39.21917 


0.19335 


2.946 


333.980 


171.178 


8.172 


51240 


(63) 


1999 


CPl33 




C92411 


39.37318 


0.27060 


9.242 


338.951 


182.862 


0.078 


51240 


(60) 


1999 


CMi58 
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Table 6 — Continued 



code 


a 


e 


i 





id 


M 


MJD 


At 


MPC name 






[AU] 




[dcg] 


[deg] 


[deg] 


[dcg] 












C98602 


41.38725 


n c\f\i r" n 

0.09152 


13.265 


10 n o o 

123.288 


n nn '^o o 

220.788 


131.912 


51600 


2 


1999 


CQl33 


i,q 


C98603 


43.76827 


n n o o o 1 

0.08381 


0.805 


n 1 1 n 

278.110 


1 n '^n A 

162.704 


36.595 


51600 


2 


1999 


CO 153 


i,q 


C98703 


42.24989 


0.0 


1.727 


1 o r" 1 o n 

135.130 


347.987 


0.000 


51220 


(6) 


1999 


CRl33 


i,q 


U99411 


44.09712 


0.04544 


2.698 


1 on n^r" 

139.275 


n r" r'no 

25.598 


onn /^ni 

329.621 


r 1 /^nn 

51600 


2 


1999 


CUi53 


i,q 


C99502 


43.36260 


n no /I n r 

0.08425 


n o /I o 

0.843 


164.676 


r" A non 

54.980 


non non 

282.282 


r 1 /^nn 

51600 


2 


1999 


CHi54 


i,q 


C04700 


42.95436 


0.0 


3.174 


ion Af\r* 

139.496 


1 O 1 o 

178.818 


0.000 


51440 


/oo\ 

(33) 


1999 


RS214 


i,q 


C04707 


42.52058 


n n /I o "1 

0.04681 


2.578 


1 o o c\rT' 

138.065 


151.541 


27.528 


51800 


2 


1999 


R'T214 


i,q 


C04910 


95.52823 


0.68200 


4.169 


1 O '7 C^f\ A 

137.804 


n 1 r" o n 

261.582 


o /I r" A A r\ 

345.449 


51460 


(68) 


1999 


RU214 




uUolUo 


45.03200 


0.07981 


1.149 


144.629 


O OO 1 

357.831 


1 on 1 nn 

180.129 


c 1 /I y< n 

51440 


(34) 


1999 


RV214 


i,q 


C09804 


43.21204 


0.07676 


1 o^n 

1.370 


n OT/^ 

0.376 


nn Ti A 

90.714 


non o 1 r" 

280.315 


r 1 onn 

51800 


o 

2 


1999 


RW214 


i,q 


C 10300 


44.90203 


n n y< 1 TO 

0.04172 


A 1 O 

4.818 


o ^no 

3.793 


c\ Ai^ nn y< 

247.094 


111 r" ^1 

111.571 


r 1 onn 

51800 


o 

2 


1999 


RX214 


i,q 


C10801 


45.55261 


n 1 o o o o 

0.18833 


13.689 


o n ^ o o o 

327.338 


^n t r A 

72.154 


o /I n rr^n 

340.779 


51800 


2 


1999 


RY214 




L.10908 


87.18745 


n CT '70 O CT 

0.57385 


20.301 




1 1 ^7 cr n 

117.659 


8.376 


51460 


(65) 


1999 


RZ214 




C13601 


AC) o "1 or\ 

48.31520 


0.34575 


7.814 


4.997 


1.548 


0.096 


51460 


(66) 


1999 


RB215 




C13604 


44.33566 


0.05693 


1 onn 

1.399 


1 on noT 

189.037 


n /I n Tnr' 

249.705 


nnr" Tnn 

295.722 


c 1 onn 

51800 


2 


1999 


RC215 


i,q 


L 14209 


120.78728 


0.68630 


25.884 


m n ooT 

210.337 


t A t n/^ y< 

141.064 


n ^n r 

2.795 


r 1 onn 

51800 


2 


1999 


RD215 




C18604 


45.16729 


0.11648 


1 A A 

1.344 


t Af\ n^o 

149.278 


1 1 n i^c\ r' 

112.725 


r' r" r' ^o 

55.573 


c 1 onn 

51800 


2 


1999 


RE215 


i,q 


C18606 


43.96413 


n n 

0.0 


3.675 


on^ r" ^n 

327.579 


0.196 


0.000 


51460 


(65) 


1999 


RF215 


i,q 


CzloOz 


47.46151 


n 1 r o /I o 

0.15843 


0.252 


294.642 


1 no o/^ o 

122.868 


ono o o /I 

298.384 


cr 1 onn 

51800 


o 

2 


1999 


RG215 




C21903 


43.79253 


0.15198 


10.212 


276.873 


68.428 


357.760 


51800 


2 


1999 


RH215 


i 


U22100 


59.78767 


n yi nm o 

0.42018 


1 n i^t n 

19.719 


OI y< noT 

314.937 


o /^no 

43.698 


o r" nn r" 

355.025 


r 1 onn 

51800 


2 


1999 


RJ215 




U 24402 


on ^nni o 

39.72913 


n 1 /^/^ A A 

0.16644 


11 A nn 

11.499 


1 OT T/^ 1 

137.761 


1 no A Af^ 

103.443 


1 1 n f\f^f* 

110.966 


r 1 /( /^n 

51460 


(65) 


1999 


RK215 




C24411 


43.21648 


n n /I nn/^ 

0.04096 


12.404 


1 /in o ^1 

140.634 


r* C) nno 

68.292 


1 n ^1 o o 

159.483 


51800 


2 


1999 


RN215 


i 


C24809 


A A m o /^ o 

44.91243 


0.10067 


1 -\ A A 

1.144 


1 o r" /^n A 

185.694 


185.986 


0.104 


51460 


(65) 


1999 


RR215 


i,q 


C24902 


^ o n o n J o 

43.08942 


0.0 


21.907 


inn o n o 

192.898 


1 'Td o A r\ 

179.849 


0.000 


51460 


(63) 


1999 


RT215 


i,q 


C25004 


43.59152 


n nn /i '7 

0.09467 


7.718 


14.298 


o 1 mo 

316.918 


36.841 


51800 


2 


1999 


RU215 


i,q 


U2ol01 


44.98636 


0.18949 


21.975 


OKI 01 o 

351.818 


27.099 


o ccn o Ac\ 
359.849 


51800 


2 


1999 


RV215 




C25103 


39.63957 


0.24625 


10.424 


253.992 


180.867 


322.404 


51460 


(63) 


1999 


RW215 




C37507 


42.05898 


0.0 


0.894 


99.349 


223.581 


0.000 


51440 


(32) 


1999 


RX215 


i,q 


C37701 


45.48866 


0.24329 


22.180 


326.622 


51.565 


328.531 


51800 


2 


1999 


RY215 




C37703 


44.04333 


0.11336 


0.563 


84.582 


297.742 


315.250 


51800 


2 


1999 


RZ253 


i,q 


C38006 


102.15245 


0.69681 


25.492 


341.669 


336.567 


0.548 


51800 


2 


1999 


RZ215 




C41409 


44.62426 


0.10334 


0.787 


192.150 


80.444 


69.830 


51800 


2 


1999 


RA2I6 


i,q 


C41607 


47.95883 


0.29806 


12.669 


175.730 


208.315 


345.437 


51800 


2 


1999 


RB2I6 




C41611 


44.40277 


0.08068 


0.588 


190.538 


344.286 


180.149 


51460 


(62) 


1999 


RC2I6 


i,q 


D17405 


43.19458 


0.18511 


4.799 


175.196 


68.292 


312.820 


51640 


(37) 


2000 


FX53 
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Table 6 — Continued 



code 


a 


e 


i 


n 




M 


MJD 


At 


MPC name 


sim*^ 




[AU] 




[deg] 


[deg] 


[deg] 


[deg] 












D30611^ 


43.62227 


0.06202 


0.565 


176.852 


301.264 


57.670 


51800 


7 


1994 GVg 


i, 


Q 


D30711 


39.40243 


0.20660 


7.164 


1.681 


243.983 


316.246 


51640 


(35) 


2000 GKi47 






D45503 


39.58196 


0.19165 


18.077 


96.123 


173.826 


292.969 


51640 


(33) 


2000 GMi47 






D45809 


41.20546 


0.00000 


28.975 


8.913 


179.957 


0.000 


51620 


(1) 


2000 GWi46 


i, 


Q 


D45904 


44.74976 


0.0 


2.718 


13.462 


134.077 


42.017 


51640 


(34) 


2000 GYi46 


h 


q 


D45906 


44.60190 


0.0 


0.664 


344.290 


204.949 


0.066 


51640 


(34) 


2000 GXi46 


h 


Q 


D49511 


43.87783 


0.00000 


1.999 


215.697 


357.395 


0.000 


51620 


(2) 


2000 GZi46 


h 


Q 



Note. — This table lists the orbits of all objects discovered, excluding the 7 lost objects which had 
insufficient timebases (2 hours) to provide meaningful orbits. The Keplerian orbital elements a, e, i, $7, 
uj, and M represent semimajor axis, eccentricity, inclination, longitude of ascending node, argument of 
perihelion, and mean anomaly, respectively. MJD is the Modified Julian Date of the orbit computation, 
and At is the timebase in oppositions (days). Orbital elements were computed independently by the Minor 
Planet Center and by David Tholen (Univ. of Hawaii) . 

^This known object was serendipitously imaged in survey fields. 

^This column identifies simulation in which the object was used: q and/or i. All CKBOs were used in 
the i simulation, and all CKBOs discovered in ecliptic fields were used in the q simulation. 



Table 7. CLF and DLF Computation 



mRrange^ iV^Lp'' ^^clf' ^ ^dlf" iVcLF^ Sdlf^ Sclf^ 



21.0 


- 21.5 


1 


1 


0.83 


1.2 


1.2 




'-'•'J -3-0.02 


03+° °^ 


21.5 


- 22.0 


2 


3 


0.83 


2.4 


3.6 


06+° °^ 


in+° °® 


22.0 


- 22.5 


4 


7 


0.83 


4.8 


8.5 




Q-i+o-ii 

'-'•^'3-0.07 


22.5 


- 23.0 


16 


23 


0.80 


19.9 


28.4 


^4+013 




23.0 


- 23.5 


18 


41 


0.70 


25.9 


54.2 






23.5 


- 24.0 


25 


66 


0.34 


73.8 


128.1 


1 QS+O-40 


o 44+0.46 
•J-^^-0.46 


24.0 


- 24.5 


7 


73 


0.11 


63.3 


191.4 


1 7n+^-'^^ 

^■''^-0.57 


r. -14+0.87 


24.5 


- 25.0 


1 


74 


0.02 


46.2 


237.5 




6.38^?:?? 



'^the apparent red magnitude range 

^the number of KBOs found within 0.5°of the echptic in the mji range 
'^the cumulative number of ecliptic KBOs found 
'^mean efficiency correction e for the given mji range 

''the bias-corrected number of KBOs, computed by summing 1/e 
(Equation ||) for all objects in the magnitude range 

^the cumulative bias-corrected number of KBOs 

^the bias-corrected surface density for the given magnitude range, equal to 
-^dlf/^) where A = 37.2 sq deg; errors are computed from 1 a Poisson errors 
for iV^LF 

'^the bias-corrected cumulative surface density, errors are summed in 
quadrature from the Sdlf errors 
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Table 8. Classical KBO Size Distribution Model Parameters 



Symbol 


Value 


Distribution 


Description 


a 


40.5 - 46 AU 


n{a)da ~ a}~^da 


semimajor axis 


pa 


2 




semimajor axis power 


e 


- 0.25 


uniform 


eccentricity 


q' 


q' > 37 AU 




perihelion distance 


i 


0-90 deg 


Gaussian, ii/2 lialf-width 


inclination distribution 


ii/2 


20 deg 




Half- Width of the i distribution 




- 360 deg 


uniform 


argument of perihelion 


n 


- 360 deg 


uniform 


longitude of the ascending node 


M 


- 360 deg 


uniform 


mean anomaly 


r 


50 - 1000 km 


n{r)dr ~ r~'^dr 


radius 


q 


fitted 




size distribution index 


PR 


0.04 




geometric red albedo 


A^CKBOs(i^ > 100 km) 


fitted 




number of CKBOs with D > 100 km 




20 




number of radius bins (log intervals) 




50 — 1000 km 




radius bin range 



"^In the circular orbit case, p corresponds to the power of the decrease in ecliptic plane surface density 
as a function of heliocentric distance R, Sed ~ R~p. 



Table 9. Selected Size Distribution Measurements of the KBOs 



rriR Range Number q Reference 

of Discovery of KBOs found 



21.1 


- 24.6 


86 


^•U_0.5 


This Work 


25.5 


- 27.2^ 


2 


3.6 ±0.1 


Chiang & Brown 1999'" 


23.8 


- 26.7 


6 


3.7 ±0.2 


Luu & Jewitt 1998'" 


23.0 


- 25.8 


5 


4-O-0.6 


Gladman et al. 1998*^ 


20.6 


- 23.0 


13 


4.0 ±0.5 


Jewitt, Luu & Trujillo 1998'' 



magnitude 

''Extrapolated from CLF slope, a, via Equation ^. 
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Table 10. Simple Inclination Model 



Ecliptic Latitude [deg] 





10^ 


20 


Number of Fields 


119 


94 


19 


Field Area [sq deg] 


37.2 


29.8 


6.0 


Number of KBOs 


74 


11 


1 


Surface Density [#/sq deg]^ 


1 on+0.22 


n s4+°-^2 

u-J^-o.io 


0.171°:?^ 


R{W°,0°) and i?(20°,0°) 




0.191^6' 


o.o9l°:J? 


H/2 [deg] 




1411 


19+f 



^Results for ecliptic latitude (3 = +10° are consistent with 
those of /? = — 10°, so were combined. 

^Error bars were computed assuming Poisson detection 
statistics (Kraft, Burrows & Nousek 1991). 
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Table 11. KBO Inclination Model Parameters 



Symbol 


Value 


Distribution 


Description 


a 


40.5 - 46 AU 


n{a)da ~ a}~^da 


semimajor axis 


pa 


2 




semimajor axis power 


e 


- 0.25 


uniform 


eccentricity 




q' > 37 AU 




perihelion distance 


i 


0-90 deg 


Gaussian, ii/2 half-width 


inclination distribution 


ii/2 


fitted 




Half- Width of the i distribution 


CO 


- 360 deg 


uniform 


argument of perihelion 




- 360 deg 


uniform 


longitude of the ascending node 


M 


- 360 deg 


uniform 


mean anomaly 


r 


50 - 1000 km 


n{r)dr ^ r~'^dr 


radius 


q 


3.7 




slope parameter 


PR 


0.04 




red albedo 


NkbOs{D > 100 km) 


fitted 




number of CKBOs with D > 100 km 




45 




number of inclination bins 




0°— 90° 




inclination bin range 



'^In the circular orbit case, p corresponds to the power of the decrease in ecliptic plane surface density 
Seel as a function of heliocentric distance R, Seci ~ R~p. 



Table 12. Minimum gout Needed 
to Explain Observed Edge 



Qout 



50.0 10 

25.0 7.4 

12.5 6.2 

6.3 5.6 



^Minimum radius [km] for which 
mass is conserved across the edge 
boundary. 



